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Abstract: Following the recent experimental realization of synthetic gauge 
potentials, Jean Dalibard addressed the question whether the adiabatic ansatz 
could be mathematically justified for a model of an atom in 2 internal states, 
shone by a quasi resonant laser beam. In this paper, we derive rigorously the 
asymptotic model guessed by the physicists, and show that this asymptotic 
analysis contains the information about the presence of vortices. Surprisingly, 
the main difficulties do not come from the nonlinear part but from the linear 
Hamiltonian. More precisely, the analysis of the nonlinear minimization problem 
and its asymptotic reduction to simpler ones, relies on an accurate partition of 
low and high frequencies (or momenta). This requires to reconsider carefully 
previous mathematical works about the adiabatic limit. Although the estimates 
are not sharp, this asymptotic analysis provides a good insight about the validity 
of the asymptotic picture, with respect to the size of the many parameters 
initially put in the complete model. 

Resume : Suite a la realisation experimentale de champs de jauge artificiels, 
Jean Dalibard a souleve la question de I'approximation adiabatique pour un 
modele d'atome a deux niveaux, eclaire par un faisceau laser resonnant. Dans 
cet article, nous derivons rigoureusement le modele asymptotique devine par les 
physiciens et montrons que cette analyse contient I'information sur la presence 
de vortex. Les difhcultcs, et c'est une surprise, ne viennent pas du termc non 
lineaire. Plus prccisement, I'analyse du probleme non lincaire, et la reduction 
asymptotique a un modele plus simple, reposent sur une separation precise 
des grandes et basses frequences (ou grands et bas moments). Cela nccessite 
de reconsidcrer avec soin les resultats mathematiques existants sur la limite 
adiabatique. Bien que les estimations ne soient pas optimales, elles fournissent 
une bonne intuition sur la validite du modele asymptotique, par rapport aux 
tallies des differents parametres initialement mis dans le modele. 
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1 Introduction 



A lot of interest, both in the mathematical and physical community, has been 
devoted in the past 10 years to the study of the rotation of a Bose Einstein 
condensate: experiments jMCWDll IMCWD2] . theoretical works (see jCooilFbt] 
for reviews), mathematical contributions |Aft[ ILiSe] . In the first experimental 
production of such rotating Bose Einstein condensates, a rotating laser beam 
was superimposed on the magnetic trap holding the atoms in order to spin up 
the condensate by creating a harmonic anisotropic rotating potential |MCWDT1 
IMCWD2| . Recently, a new experimental device has emerged which consists in 
realizing artificial or synthetic gauge magnetic forces and leads to the formation 
of vortex lattices at rest in the lab frame jLCGPS) . A colloquium jPGJOj has 
analyzed in detail the artificial gauge fields and their manifestations. In order 
to understand the main ingredients of the physics of geometrical gauge fields, 
|DGJOj (see also [GCYRD] ) study the case of a single quantum particle state 
with a 2 levels internal structure. More complex systems with more than 2 
internal levels are also discussed in jPGJOj but we stick here to the simpler 
case of 2 levels, which contains all the mathematical difficulties. A key issue 
is to determine whether one internal state can be followed adiabatically. A 
question raised by Jean Dalibard is to analyze in particular whether vortex 
formation may break down the adiabatic process. In jPGJOj . some conditions 
are provided, that we want to analyze from a mathematical point of view. 
We are interested in the minimization of the energy 

= I |V0|2 + K(a:,y)|0P + ||0|''da;dy 

+^i,,,^^x){cp , -cos(0,J.t)) 
where = (tt'^J^ ^ ^nd \Hx,yr = \M^,y)\' + I02(x,2/)p. We 



h{x,y)^ 

prescribe that the norm of is 1. Here 0i and (j>2 are the internal degree 
of freedom of a particle: ground and excited state of the atom. It is assumed 
that the atom is shone by a quasi resonant laser beam. The functions Q^^g^^x), 
(fik (y) and 9g^ {x) are given as in |DGJOj by 



n^,,e^{x) = K^{ — ) with n{x) = \n + x^ , 

Oi. (x) = ^(f ) with cos(^(x)) = , sin(^(a;)) = , 

fk{v) = ^{ky) with Lp{y) = y , 

where is the phase of the propagating laser beam while 9_ is the mixing angle. 
We define 

,w ^ r,, .( cos{e_{x)) e^^(«) sin(0(x))\ ,^ 
sm[U(x)) ~ cos[Oxx)) J 
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The matrix M models the couphng between the atom and the laser. The 2x2 
matrix M{^^=, y/ljty) can be diagonalized in the bases {ip+^tp-) respectively 
associated with the eigenvalues ^^{x), 

with C = cos [\e{-^)^, S = sin (^i^(-^)j and ip = -/I^y. When the 
particle follows adiabatically the eigenstate ij}- > this corresponds to set formally 
(j) = u{x,y)ip-, where u{x,y) € C in f^. Then u minimizes a Gross-Pitacvskii 
type energy functional with a modified trapping potential called the geometrical 
gauge potential. The scalar potential VK,{x,y) = VK,{x,y)ldc2 will be adjusted 
in order to produce a harmonic potential after the addition of the geometrical 
gauge potential from the adiabatic theory. We want to justify the adiabatic 
approximation for states close to ip_ and analyze the error term between the 
initial and effective Hamiltonians. 

After a rescaling, the parameter occurring in the experiments have the fol- 
lowing orders of magnitude 

K - 10^ , G - 600 , 4-25 , /c - 50 , 

but other values can be discussed. Conditions on the strength and spatial extent 
of the artificial potential have to be prescribed in order to induce large circula- 
tion. Two cases, i^k > 1 and i^k < 1, can be distinguished and the problem 
has to be rewritten in two different ways in order to apply semiclassical tech- 
niques. In fact, we will focus on the case (^k > 1, corresponding to the previous 
numerical values. The complete analysis is carried out in the asymptotic regime 
£K.k — > -l-cxD but some partial results are also valid for (^.k < 1 or . 

A change of scale (j){x, y) = ^^"0(^/^2;, \J^y) yields a new expression for the 
energy: 

+ ^^(^)(^, { - rj^^l*^^ e-(^^)sin(,(^))N 

^Vfc4^^'^'l^e-^^(V^'^)sin(^(^)) _cos(^(^)) j"^'^ 

+ ^Ji;\Uxdy. 

According to the two cases (^.k > 1 or £„fc < 1, we define a small parameter 
e that allows to rescale the energy. In fact, we define rather the parameter 
e^+^'^j where 6 can be taken as a first step equal to 5/2. The exponent S > 
is a technical trick which provides the right quantitative estimates for the 
adiabatic approximation with a quadratic kinetic energy term . The suitable 
choice of this new parameter S is discussed further down in this introduction, in 
Subsections ll.ll and ll.2l The small parameter e > is thus introduced according 
to the two cases: 
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if > 1, then 



1^ l^f^K. f^^K. 



(This leads in our example to e2+2'5 = 2.5 10 '^.) 
if < 1, then 

e^+^' = ^,5>Q , G, = ^ = Gk£,s'+^'. (1.4) 

We define 

r-ir r)-{ ^^'^^ if > 1 ' n 5) 

■'-^■'^^■'y)-\ iffc^«<l. ^^-^^ 

In both cases, this leads to 

«:-if„(0) =£-,(^) (1.6) 



where 



and 



- {i' , Huni^) + %^ y IV'I' dxdy , (1.7) 



i/L»» = -e"V,Ty£^A + y,,,(x,y) + A/(,/-x,./^2/), (1.9) 



ye,,(x,y) = ^-W^{J-^x,J^y) (to be fixed). (1.10) 
The quadratic energy (linear Hamiltonian) is 

£<iuad,eW = / e'+'V,rjVV|' + K,r(a;,y)|V'|' 

+ , Mi l^x, .f^y)i^)c2 dxdy (1.11) 

= {i,,HLrn^). (1.12) 

At least when = Tj, = 1 and (5 = 0, this problem looks like the standard 
problem of spatial adiabatic approximation studied in jSor] . |MaSoj . jMaSo2j 
|PSTj ■ although it requires some adaptations because the symbols are neither 
bounded nor elliptic. 

We shall consider the asymptotic analysis as £ — > with uniform control 
with respect to the parameters G, which allow to fix the range of validity 
of the reduced models. Then we shall consider the asymptotic behaviour of 
the reduced model and the whole system large. Specifying the right 

assumptions on Ve^rix^ y) or V„, possibly in a scale depending on (^k, fc), is also 
an issue. 
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1.1 Main result for the Gross-Pitaevskii energy 

We shall choose the potential V^^t in (|1.10p such that after the addition of the 
adiabatic potential in the lower energy band, the effective potential is almost 
harmonic. Namely, we assume 



^2+25 



1 



(1.13) 



with the potential v chosen such that 

«(x, y) = (x^ + y^)xv{x^ + y^) + (1 - Xvix^ + 2/')) (1.14) 
X. eCo°°([0,2)) , 0<x. <1, X. = lon[0,l], 

and £v > parametrizes the shape of the quadratic potential around the origin. 
With these assumptions on the potential V^^t, if one chooses ip = u{x,y)ip^^ 
where ip- is the eigenfunction corresponding to —rt{xy^ ^xl^y) of the matrix 
^{x^TxlTy)M(x^ Tx/Ty, yyjTy/Tx), then the linear part Hun in £e is formally 
replaced by the scalar e^'^^^TxH- where 

H- = -dl - (dy - z— =^=^) + -2^w(^/^x, ^y) . (1.15) 

\ Zy 1 + TxX / t^yix 

In the limit r^; — >■ 0, a natural (e, T)-independcnt scalar reduced model emerges: 



£h{i 



-dl {dy y )^ + 



G 
2" 



u{x,y) e C. 



Because k is large, or e is small, it is natural to expect that the ground state of 
£s is close, up to a unitary transform, to a vector u(x,y)ip- . This is the aim of 
the adiabatic theory and leads to a scalar problem. In order to get good bounds 
on the energy, we need to study the limit Tx small at the same time. 

In all our work £v > and G > are assumed to be fixed, while the 
asymptotic behaviour is studied as e — >■ and — > . 
The quadratic part of the above energy is associated with the Hamiltonian 

H,, = ^di-idy-'-^r+^^±y^, {o<ev<+oo), (1.16) 

with the domain 

n2^ lueL^iR^), J2 Ik^^^^^lU^ < +^ > , q={x,y)eR\ (1.17) 

[ l"l + l/9|<2 J 

endowed with the norm H/H^^ ~ E|q|+|/3|<2 Ik"^? '"IIl^ ^'^"-^ 1^'-' corresponding 
distance d-H2 ■ It is not difficult (see Section |4?2|) to check that the minimization 
of £h{'p) under the constraint H'/'IIl^ = 1, admits solutions and that the set, 
Argmin Eh , of ground states for £h is a bounded set of '^2 . 



5 



Definition 1.1. For a functional £ defined on a Hilhert space % (with +00 as 
a possible value), we set 

£min = inf £{u) and Argmin f = {u ^H, £{u) = £min and\\u\\L2 = 1} . 

Il'"llt2=l 

Theorem 1.2. Fix the constant £v,G and S and assume, for some C's = 

5 

e^'<^. (1.18) 

Let X = (XI7X2) be a pair of cut-off functions such that Xi + X2 ^ ^ '^^ ) 
Xi G C^(M^) and Xi = 1 i'^ a neighborhood of . 

There exists vq = i^ev,G S (Oi ^] o,n-d for any given S > 0, there are constants 
Ts = t{£v,G,S) > 0, Cy^,s = C{£v,G,6,x) > 0, and a unitary operator U = 
U{e,T,£v,G,d) which guarantee the following properties 

• The energy £^ introduced in (|1.7p admits ground states as soon as < Tg, 
with 

If . _ p2+25 c . I < 2+25 i 

/ - " \ 

^ I e ^v^ ^ Q I \ 

• For any ip € Argmin £^ written in the form ijj = U \ v , the 



vector a = [ ] satisfies 



\a+\\L2 <Cse^+^'T, , \\a\\H'^<^ , + M^e < Cs . 



||X2(T-f'.)a-||L2 < Cx,5Ta? , 

1 2:.o 

c^«2(xi('r^'-)a-, Argmin £h) < C^^s{tx^ + e) : 



lla+lk- < Cse^+^ , dL=(xi(Ti.)a-, Argmin £h) < Cx,5(tJ^" + e)^ . 

All the constants can be chosen uniformly with respect to 5 ^ (0, Sq] for any 
fixed i5o > . 

Remark 1.3. The proof is made in two steps: 1) the limit e — > corresponds 
to the adiabatic limit for the linear problem and allows to replace the linear part 
Hhint of£e, by the scalar e^'^^^T^H- given by p.lSp ; 2) the limit ^ allows 
to reduce the asymptotic minimization problem to a simpler one where the linear 
Hamiltonian is exactly H^^, given by (|1.16p . 

One main point in the proof is to have precise energy estimates for the limiting 
problem. In our case, we obtain them in the limit — )■ 0, because explicit 
calculations are more easily accessible when the magnetic field is constant, that 
is in the case of £H.min ■ In theory, if one had precise energy estimates for a 
general for the intermediate adiabatic model with the linear part H- , complete 
results could be performed for general values of . 



6 



Remark 1.4. The exponent vq = 1^1^,0 is a Lojasiewicz- Simon exponent (see 
Remark \4-3\ for an explanation and a. e. fLoj] \BCRf for the definition of Lo- 
jasiewicz exponents and \Simf for its extension to PDE problems). It is i when 
the Lagrange multiplier associated with u G Argmin £}j is a simple eigenvalue 
of + G\ijj\'^ . There are reasons to think that it is the case for generic values 
{£v,G) € (0, +00) X [0, +00), namely outside a subanalytic subset of dimension 
smaller than 1 (and possibly 0). Nevertheless for G ~ 0, there is a discrete set 
of values of iy for which He^ has multiple eigenvalues. 

Consequences and applications: One issue is whether the presence of vor- 
tices (zeroes of the wave function with circulation around them) might break the 
adiabatic approach. The answer contained in our theorem is that vortices of iJj 
in the original problem and vortices of the ground state of the Gross-Pitaevskii 
energy £h arc close and that the minimization of £h provides all the informa- 
tion on the defects of i/;. Indeed, the smallness of a+ in L°° indicates that the 
vortices of ip a-nd a_ are close, and the last estimate of the theorem provides 
that the vortices of a_ are close to that of the Gross-Pitaevskii problem. 

Numerically, in |GCYRE)] . the authors observe vortices in a system with 
artificial gauge as presented in |DGJO| and modeled with the energy They 
check that the vortex pattern is close to that of the Gross-Pitaevskii energy. 
If one wants to use our results, one may process in the following way: once 
G is fixed, choose iy such that the minimizers of £h have vortices (detailed 

conditions will be given in section IT2|) . Then take Tx > and e > small 

1 

enough so that the L°° norm of a+ and the L°° distance from x(ra?.)a_ to a 
ground state of £h arc small. In the above result, the constants are not explicitly 
controlled and this control is worse and worse when S increases. So it is not 
explicit for given numerical values of the parameters £v,G,Tx,e or equivalently 
£vj G, K, k, ■ Nevertheless it provides a framework for numerical simulations, 
where the observation of vortices can be confirmed by decreasing e and r^. . The 
parameters of the experiments are just at the border where our constants may 
become too large to provide a reasonable approximation. 

The definitions (|1.3p of e and (|1.5|) of , transform the condition p.l8|) into 

The larger 5, the better, but S must not be too large because of the value of 
the constants Cs, Cs^x- ^ value of a few units for 6 does not affect them too 
much. Another reason for keeping S small is that the initial small parame- 
ter is e^+^'^, while the final error estimate of d^^^ (xi('''i' .)a_, Argmin £h) (or 

1 

dioo (xi(tJ' .)a_, Argmin £h)) contains also an 0(e) term. 
As an example for 5 = |, the above relation becomes 
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A given precision of order e + Tx , is more easily achieved by taking k small and 
= 7^ large (for example k = 0.1, £k = 500 x 25 is better than the values 
k = 50, £k = 25 given in the introduction). Note also that the external potential 
Ve^T, defined in (|1.13p must be adjusted up to the order £2+25 _ k_ ^ 



1.2 Gist of the analysis 



Following the general idea of the founding articles |BoFo[ |BoOp| of Born. Fock 
and Oppcnheimer, it is well known in the physics literature that a Hamiltonian 
system 



{eDgf+uo{q) 



E+iq) 
E^iq) 



is unitarily equivalent to a diagonal Hamiltonian plus a remainder term; 



iD,-Aiq)r 








\Xiq)\' 



Rie) 



where {Dq^^Ai{q)) are the covariant derivatives {^A{q) is the adiabatic connec- 



tion) associated with the fiber bundles UQ{q) 



and Uo{q) 



and |X(g)|^ 



is the Born-Huang potential. 

Since Hit] and until recently ( [NeSol IMaSol IPSTI IPST2| ). this problem has 
been widely studied by mathematicians, and the remainder term is formally: 

Rie)^^s^a,iq)isD,^)ieD,^) + Oie'). 

It is smaller than for bounded frequencies (or momenta) but it has the same 
size as the main term for a typical frequency of order 7 . For a nonlinear problem 
or without any information about the frequency localization of the quantum 
states, it is important to estimate the error terms in the low and high frequency 
regimes. 

Introducing S > allows to obtain at the formal level 



.2+2<5 



{D,-Aiq)r 






iD,+Aiq)y 



1^(9)1' 



.2+4(5 



a,{q){eDg^)ieDg^) + Oie' 



3+45 \ 



where the remainder term is now 0{e^'^'^^) at the typical frequency j- and thus 
0{e'^^) times the size of the main term. Such an error estimate can be made 
in the L^-sense when applied to some wave function ip lying in {\p\ < r}, with 
p quantized into eDg, or more precisely fulfilling ip = x(eDq)ip for some x ^ 
C^{{\P\ < 2r}) with 

\\R{em < C^ie'^^' + CAr£^)|| VII < iC'^^se'')e'+''\m • 

where Cat essentially depends on the estimates of A^-derivatives of UQ{q), E±{q) . 
For a fixed (5, we choose N >2 + A5 . 
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Remark 1.5. About the choice S > 0, another strategy could be considered 
in order to optimize the exponent d, w.r.t e: under analyticity assumptions or 
m^ore generally assumptions which lead to an explicit control of Cn in terms 
of N, one could think of optimizing first Cn£^ w.r.t to N according to the 
methods of fNeSoi ISori \MaSof . As an example with Cn < Nl, this would lead 
to Cn£^ < Ce~~ after choosing N = N{e) = [^], with some C < 1 . Then 
taking 5 = S{e) ~ ^JTiagie) '^^''^^^ lead to 

CCx(e2+4«(^) +e"*) < 2CC^e-^\i^\\ . 

When (5(e) — >■ oo and as compared with e2+2i5(e) considered as the order 1 term, 
an remainder term is almost of order 2 . 

We do not consider this optimization of 6 w.r.t e, with lim£_j.o = +oo, be- 
cause our initial small parameter is e^"*"^* and the final result of Theorem \1.2l 

(the estimate about d-u^ (xi ('He 0^- 1 Argmin Eh) in the nonlinear problem), con- 
tains an ©(e) . Hence we keep 5 > Q independent of e. This is why S is still 
present in the constants of Theorem \l.S[ 



Wc need to perform a frequency (or momentum) truncation. We decompose 
a general ip into x{sDq)^p + (1 — xi^Dq))^ and use rough estimates for the part 
(1 — xi£Dq))ijj which wiU be compensated in the minimization problem for £^ by 
good a priori estimates for the norm || {l — x{sDq))tp\\ of the high-frequency part. 
We also have to check that the unitary transform U, implementing the adiabatic 
approximation, does not perturb too much the nonlinear part of £e{i^) . 

In our case, the limit e — > leads to the Born-Oppcnhcimcr Hamiltonians 

which, in a second step, in the limit Tx — > 0, leads to H^^, (with the sign —ix 
for the lower energy band). This means that the convergence to the Born- 
Oppenhcimer Hamiltonian as e — >■ has to be uniform w.r.t G (0, 1] . This 
last point requires to reconsider carefully the work of jPSTj by following the uni- 
formity w.r.t T of the estimates given by Weyl-Hormander calculus ( jHorllBoLej ) 
for T-dependent metrics which have uniform structural constants. 
This is done for the low frequency part in Section [2] while the basic tools of 
semiclassical calculus are reviewed and adapted in the Appendix |Al 
In Section |3] the error associated in the high-frequency part is considered, as 
well as the effect of the unitary adiabatic transformation on the nonlinear term. 
Once the adiabatic approximation is well justified in this rather involved frame- 
work, the accurate analysis, as well as the comparison when Tx is small, of the 
two reduced models (the one with Hi^, and the one with _ff_) is carried out 
in Section 21 This follows the general scheme of comparison of minimization 
problems: 1) Write energy estimates; 2) Use bootstrap arguments and possibly 
Lojasiewicz-Simon inequalities in order to compare the minimizers in the energy 
space; 3) Use the Euler-Lagrange equations in order to get a better comparison 
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in higher regularity spaces. 

In Section [SJ aU the information of the previous sections is gathered in order to 
prove Theorem 11.21 : existence of a minimizer for Eg in Proposition [531 key en- 
ergy estimates in Proposition 15 .41 and bounds for minimizers in Proposition [521 
Some comments and additional results are pointed out in SectionlHl namely: 1) 
the question of the smallness condition of e w.r.t appearing in Theorem 11.21 
2) the possible extension to anisotropic nonlincarities (one would have to check 
that the unitary transform implementing the adiabatic approximation does not 
perturb the nonlinear part); 3) the minimization problem for excited states, i.e. 
locally and approximately carried by '0+ instead of ip- ; 4) the extension to the 
problem of the time nonlinear dynamics of adiabatically prepared states. 

2 Adiabatic approximation for the linear prob- 
lem 

In |PST| . the adiabatic approximation is completely justified for bounded sym- 
bols or when global elliptic properties of the complete matricial symbol allow 
to reduce to this case after spectral truncation. Unfortunately, it is not the 
case here, because the eigenvalues of the symbol of the linear part HLin are 
e^^T-a:Ty|pp -1-14,1- (a;, y)^^{^j^x) ■ In |Sor| . the adiabatic theory for unbounded 
symbols is developed after stopping the complete asymptotic expansion in an 
optimal way, under some analyticity assumptions, but this would be particu- 
larly tricky here with divergences occurring both in the momentum and position 
directions. We shall see that the sublinear divergence in position makes no dif- 
ficulty after using the right Wcyl-Hormander class. The quadratic divergence 
in momentum, with the kinetic energy TxTy\p\'^ is solved by first considering 
truncated kinetic energies and using the additional scaling factor e^^ , 6 > 0, in 
front of the kinetic energy term. 

Our problem shows an anisotropy in the position variables (x, y) . The analysis 
of the linear problem can be treated in 'R'^ . Then wc split the position and 
momentum variables, q eW^ and p gW^, into: 

q={q\q") , p^b',^") q' , p' E R"' , q' \ p" E R"" , d' + d" = d , 

and the pair {tx, Ty) is accordingly denoted by (r', r") e (0, 1]^ . 
From this section, some notions and notations related with semiclassical anal- 
ysis arc used. In particular, the notation Suirrir, gr) refers to classes of {e,t)- 
dcpendcnt symbols of which the seminorms are uniformly controlled w.r.t to the 
parameters (e, r) G (0, Eq) x (0, 1]^ . For accurate definitions, we refer the reader 
to appendix [^ where all the necessary material is reviewed and adapted for our 
analysis, assuming knowledges about Frechet spaces and generalized functions. 
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2.1 Born-Oppenheimer Hamiltonian 

ailtonian in = H{ 
H{q,p,e) = £2Vr"|pp^(r'r"bn + V{q,T,e) 



Consider the Hamiltonian in = H{q,eDq,e) with the symbol on M^'^ 



e''T'r"\p\'j{T'T"\p\')+uo{q,T,e) 



E+{q,T,e) 
E_{q,T,e: 



ul{q,T,e) , 



with i5 > and (t', t") £ (0, 1]^ . The following properties, with the splitting of 
variables q — {q',q") € K"*, are assumed: 



E± e Si, 



E+ {q, r, e) - E_ {q, re) > C^^^l —q^ ) , 

t" 



(2.1) 



V9"';-^2(c)) 



7 e C(^(R; M+) with 7 = 1 e [0, 2r^^\ . 

With these assumptions we are able to justify the Born-Oppenheimer adiabatic 
approximation for (5 > 0, t',t" G (0, 1]. We shall work with the r-depcndent 
metric 



5r 



on the phase-space R ^^p, w hich is checked to have uniform properties w.r.t t S 
(0, 1]^ in Proposition I A. 91 . The exact definition of the parameter dependent 
Hormander symbol classes, Su(rn, gr', M2{'C)), is given in Appendix |^ (see in 
particular its meaning for r-dependent metrics in Appendix IA.4[) . We shall use 
the notation 

1 



'^q'){./7^pY 

with the meaning of the exponent oo being the same as in C°° of 0{e°° ) . 
Theorem 2.1. There exists a unitary operator U = U{q,eDq,T,e) with symbol 
U{q,P,T,e) = uq {q, T, e) + eui (q, p, r, e) -t- e^us {q, p,T,e), 



{JTiTq'){V^py' 



such that 
U*HU = 



hBO.+ {q,£Dq,e) 








hBo-{q,£Dq,£) 
e^+^^Riiq, eDq, r, e) + e^+^^R2{q, eDq,T, e) , (2.2) 
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with hBo±{<l,P,''', s) equal to 

= £^^T'T"i\pTeA\^ + \eX\^)+E± 



.k=l 



(2.3) 



when \Jt't"\p\ < r^, and 



Xk —Ak 



The remainder terms satisfy Ri, R2 € S'u(l, (jr; A^2(C)) and R2 vanishes in 
and those estimates are uniform w.r.t S G (0,So\ (and r S 

(0,1]2, £G (0,£o];. 

Remark 2.2. • The remainder term is really negligible only for 6 > 0. This 
is explained in Subsection \2.4\ 

• The (Weyl)- quantization of e"^^ t' t" {\p^ £A\'^ + \£X\'^) + E± is nothing but 

d 



Y.^dg'-TiAkf + \Xu 



.k=l 



E. 



± ■ 



2.2 Second order computations for space adiabatic ap- 
proximate projections of the reduced Hamiltonian 

We shall consider the matricial symbol, on M^'^, 

Hiq,p, T,e) = f,{p, t) + E+{q,T,£)nn{q,T,£) + E^{q, T,e)(l - no((7,r,e)) 

where ^oil, t, e) ~ tto{q, t, e)* ~ ^oil, t, e)^ e M2{C) , 7, E± real-valued, with 
6 > and the following properties: 



E±eSM—q'),g,r) , Ilo€ Suil,g,.r;M2iC)), (2.4) 

V T 

^dq'^ t" ,,2 
with gq^r = I— 1- —dq , 



E+ [q, T,£)~E- {q, T, e)>C-\J ^q') , 



(2.5) 



f,{p,T)^£''fi{VT^p) , fieC^{R'';R) , S>0. 



For conciseness, the arguments p, q and the parameters r, £, will often be omitted 
in d^f.ip) and d^E±{q) or d^Uo . 
According to Appendix I A. 41 the metric 

5r = 59,r + , = ^= + -fdq ' 



{sJV^p)' 
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has the gain function 



X{q,p) = mm — ^ , > {\Jt't"v) 



The e-quantized version of the symbol A^ici^-p) wiU be denoted 

Note that the symbols -E+ — E- is elliptic in its class Su{{\J~^q') ■, Qt) , and 

1 



{Vt't"p)°° 



Our aim is to compute accurately the adiabatic projection II'") {q, p) = no((7, £) + 
eni(g,p, £) + ••• + £"n„(g,p, e) such that 

n(") ofK") = nH +c'(£"+i) ^ r7?,n(")j =c'(e"+i). 

The general theory presented in [PST| . tells us that the asymptotic expansion 
can be pushed up to n = oo, but we will do here accurate calculations up to 
n = 2 (with additional information for n = 3) and then discuss the influence of 
the factor e^^ . Those are feasible and rather easy because the kinetic energy 
term and the two-level potential are simple. This allows to reconsider accurately 
the arguments sketched in [PSTj for the Born-Oppenheimer case with all the 
technical new peculiarities of our example. Note also that in |MaSo2| chapter 10 
explicit calculations have been made up to order n = 5 but in the slightly dif- 
ferent framework of time-dependent Born-Oppenheimer approximation oriented 
to polyatomic molecules : no use of the exponent (5 > 0, no divergence as g — > oo 
and no ellipticity problem, no extra-parameter r and the techniques are slightly 
different although still relying on semiclassical calculus. 

Like in |PSTj . (Sot], |MaSo| . |MaSo2] . the calculations are first done at the 
symbolic level and we write 

{A[e) - Os(£'^)) ^ {e-^A G S,(l, g,; X2(C))) . 

For a matricial symbol A{ci^p) (possibly depending on (£,t)), it is convenient 
to introduce the diagonal and off-diagonal parts 

^^(g,p) no(g)A(q,p)no(g) + (l-no(g))A(q,p)(l-no(9)) (2.6) 

AO^(g,p) no(<z)A((z,p)(l - no((?)) + (1 - no(g))A(q,p)no(g) , (2.7) 

where we recall no(g) = no((7,T, £) . Note the equalities 

(Ai?)^ = A^i?^ + A^^ijo^ , (yli3)0^ =A^50^+A0^5^. (2.8) 
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The "Pauli matrix" 

(73((?,T,£)=2no(g,T,£)-Idc2, (2.9) 

will also be used, with the relations 

al {q) = Idc2 , a3{q)A^{q,p)a3{q)^AD{q,p), 
asiq)A^''iq,p)asiq) = -A^''{q,p), 

a3{q)AO°{q,p) = Uoiq)Aiq,p){l - Roiq)) - (1 - Ilo{q))Aiq, p)noiq) . 
We are looking for 



n(") iq,p, T, £) = ^ em,{q,p, T, e) € 5„(1, gr; A^2(C)) 

3=0 

with n, e 5„(1,5,;A^2(C)), 



and such that 



(2.10) 
(2.11) 
(2.12) 



Like in |MaSo| . |PST| . this system is solved by induction by starting from 
n^") (q, p, T, e) = Ho ((?, T, e) , with 



Gn+l 



.-("+!) 



n) 



mod Cs(£) , 



(2.13) 
(2.14) 



F, 



n+l 



-(,1+1) 



-(n+l) 



iff (n(") +e"+in,f+j 

-(n^") + £"+in^+i)f if] mod Os(e) 



n: 



OD 
n+l 



E+{q)-E^{q) 



CT3^„+1 - 



mod Osis) , 



1 



E+{q)-E^{q) 



(2.15) 
(2.16) 



The general theory says that the principal symbol of F„_|_i is off-diagonal, 
Fn+i — Ff^^i mod Os{£) , and F^+i can be chosen so that 



(2.17) 



Below are the computations up to n = 2 in our specific case. In these com- 
putations, we shall use Einstein's summation rule Skt^ = s/jt*^ with the 



coordinates {pk,q ) or {p ,q ) with p*^ 



e.k ~ Pk like in the examples 



d 



\p\' = p'pk = PkPiS'-' = p''p'Sk,i , (9p/e).a, = {dp.f,)S'^-'— = {dpje)dg. . 
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n = 0: Start with n(°) = Ua{q, r, e) and notice dpUa = and dqUo = {dqUa)°^ . 
n = 1: Take 

Gi = Ho o Ho - Ho = , 
and nf = . 

Next compute 

[Hfnoiq, e) - no(g, e)fH\ = e"! [Af no(g, s) - no(g, e)fh] 

= --idpJ^dgkUo mod 05(e), 

and take 

and n(i) = Ho + -£!^^a3(9,.no)°^ , 

with 

n = 2: Consider now 

£2^2 = n(i)f n(i) - n(i) mod ©^(e^) 

= nofno-no + £(nofni + nifno)-ni + £2nifni mod ©^(e^) 

According to (jA.3|) . with Ilottllo ~ Hq and with 

HoHi + Hi Ho = <n?^ + nf ^n^ = Hi , 

we can take 

G2 = - [-9,.no<9p,ni + dp,Uidg.Uo] + uj . 

The first term —d„kIlodpi_Ili, with Einstein's summation rule, equals 



id^ f 



-0-3(9 fe Ho) (a f Ho) 



while the second term +{dp^Ili){dqkIlo) gives the same result. 
The third term is given by 

tt2 {(^Pkfe){dpje) ^ ^ \OD^ ( tt ^OD 

i^i = (£'_^ - £• )2 ^^(c'g'-Tioj cr3(d,aio) 

~ (^+ - £;_)2 (^g''^^o)(Cg^iloj ■ 
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Hence the diagonal second order correction is given by 



{E+ - E^fn^ = -{E+ - E^fa^G2 

= - [{dpj,){dpj,) + {E+ - E^)(dl^J,)as] (9,. Ho) (a,. Ho) 



and satisfies 



e-^'n^ e Su(^^ — ,9r;M2{C) 

Consider now 112'^: By referring to (j2.17p . F2 can be chosen as the off-diagonal 
part of [ffji^n'^^^ — n'^^^tJ^iJ] which, according to the previous steps and 
(|A.3p . equals 



1 

8 

1 



+ — [dp^^HdqkHi - dqkHdp^Ui - dp^UidgkH + dgkUidp^H] mod Os{e) . 

(2.18) 

Since dp^ p^H = {dp^ pje) as a scalar symbol commutes with d^k^J^o , the first 
term of (|2.18p vanishes. 

Similarly the factor dqk Hi appearing in the second term of (|2.18p contains three 
terms 



dqk{E+-E^) u^^AMlL 

i{dpje) 



{E+-E-) 

where the second one is diagonal (Remember (73(5, T,e) = 2no((7,T, e) — Idea) . 
Since dp^H ~ dp^fe is diagonal, we get 

j^[dp,Hd,kn^ + d,kn,dp,H]'"' 

= idpje)idpj,)a, {d,k [iE+ - E^)-'d,eU„]f'' . 



In the quantity —dqkHdp^IVi — dp^WidqkH , the derivatives 

Opk^h = (j;^ _ E ) ^3(ggtllo) 

are off-diagonal factors, while 

dqkH = {dgkE+)no + {dqkE^){l - Ho) + {E+ - £;_)(9,.no)°^ . 
With the two equalities, 

(aqfcS+)nocr3(a,/no) + cr3(a,fno)((9,fe£:+)no = {dqkE+)aj,{dqtTiQ) , 
{dqkE^){i - no)c73(a,.no) + c73(a,*no)(9,.S-)(i - no) = {dqkE^)a^{dq,no) 
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we get 

This leads to 



f 



-idg.{E++E_))a:i{dgmo) 



{E+-E^] 



and 



n: 



with 



OD 



iE+-E_) 



{dg. [{E+-E^)-'dgmo]) 



OD 



f 



2{E+-E^Y 
1 



(a,.(£;+ + s_))(9g.no), 



We have ahuost proved the 

Proposition 2.3. The pseudodijferential operator Tl{e) ~ Il{q,eDq,e) given by 
n(9,p,e) = Uo{q,e) + £ni(g,p, e) + e^n2(g,p, e) 

\OD _ ^(^Pkfs 



with Iii{q,p,e) =\li{q,p,ef {^e+ - E ) 
n2 (g, p, e) = n2 (g, p, e)^ + (g, p, e)^^ 



CT3(<9„icno) , 



n2(g,p,e) 



[(5p./e)(9pje) 



(a,. [{E+-E^)-^d,.Wo\) 



OD 



and e-'^^ni,e-^^n2e Su 
satisfies 



-,5r;A^2(C) 



{J^q'){V^pY 



Hon = n + c'(e 



3+2(5 \ 



n* = n 



3+2(5\ 



in £(L^(M''; C^)) . Moreover the estimates in £(L^(]R''; C^)) of the remainder 
terms do not depend on the parameter t = {t',t") G (0, 1]^ and S > 0, as soon 
as e € (0, £0) . 
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Proof: The above construction gives immediately 



noii = n + o{e^) , n* = n 



The first improved estimates come from the fact that 

noH-n 

contains only terms which are Moyal products with a IIi or a 112 factor, with 
cancellations up to the £^ coefficient. Both of them have seminorms of order 

For the last one, this is a similar argument after decomposing 

□ 



The above result can be improved after considering what happens at step 
n = 3 when f^^rip) = fi{\/ t't"p) is at most quadratic w.r.t p in some region. 
Before this, let us examine the remainders of order 3 in e. 
n = 3 : The remainder term 

e^Gg = n(2)fn(2) -n(2) 

= e (Hof Hi + Hif Ho - Hi) + (Hif Hi + Hof Ha + Haf Hq - Ha) 
+e3(nif n2 + Haf Hi) + e^Ii^n . 



Using the construction of Hi and 1X2 and the fact that e ^'^IIi and e ■^''112 be 



-25t 



long to 



- (y 777'?'>(vt't"p) = 

(|X2|) -(|X3 l) tells us 



5r; A^2(C) ) , the expansion of the Moyal product 



n(2)jsn(2) _n(2) = [AG(9,'no,9|ni,e) + i?G (5,00,9^02, e)] +e^+^'RG 

(2.19) 

where Rq G 5,4 — 1= — ^ ,5r;A^2(C)) and AG{.,e), Bg{-,s) have an 

V (^^a')(VT'T"v)°° I 



asymptotic expansion in terms of er, of which all the terms are bilinear differential 
expressions of their arguments. For the commutator with H, write 



iff n(2) - n(2)f ij 



= [A(p)fno(g)-no(g)f/ap)] 
{E+i\o + E^{i- no))f n(2) - n(2)f (i;+no + e^{i- Hq)) 



After eliminating all the terms which are cancelled while constructing Hi and 
112, the contributions of all three terms of the right-hand side can be analyzed. 
The contribution of the third term is similar to what we got for G3: 



£3 [AH{dl{E±IVo),dlni,e) + BH{d,{E±I{o),dpIV2,e)\ +e^+^'RH^ 
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with Rh.1 e S,(^^=i^,.gr;A^2(C)) 



With the uniform estimate of e ^^Hi, 



e ^*n2 and e "^^ fe, the contribution of the second term is estimated as 



The contribution of the first 



with Rh,2 e Su( — r=, \ 

term is £'^C//(9p/e, S^IIo, e) and we get 

AHid^giE±no),dlnue) + BH{dgiE±no),dpIl2,e)] +e^+^^RH (2.20) 

where the expansions of Ah, Bh and Ch w.r.t to e have terms which are bihnear 
differential expression of their arguments, and the remainder Rh belongs to 



Proposition 2.4. With f^{p) £^'^7i(VtV'p), assume that the third differ- 
ential dpfi vanishes in {\p\ < r} and fix r' £ (0,r) . Then the remainders of 
Proposition \2.3\ equal 



n o n - n = e^+^'i^ci + e^+^^RG..2 



-3+25 

RH,i+e'+^'RH. 



(2.21) 
(2.22) 



where RG,ior2 belong to OpS. 



1 



gr;M2{C)), RH,ior2 belong 



to OpSu 



r' r" p) " 



,5r;A^2(C)) and RgotHa = m ||Vr'r"p| < r'j . Thost 



the symbol 112 is the sum of a linear expression of and quadratic expression 



estimates are uniform for r G (0, 1]^, (5 G [0, (5o] and e G (0, £o) 

Proof: After noticing that the symbol Hi is a linear expression in dpfe while 

off'^, 

in 5p/e, the identities (|2J9)) and ((2?20)) imply 



for an arbitrary large G N, in < r| . ChooseS N > 25 and 

take a cut-off function x G ^0° (ibl such that x = 1 in a neighborhood 
{\p\ ^ f'} ■ Writing for the symbol 

s = n(2)f n(2) - n(2) or 5 = [i?f n^^) - n(2)f h 



S* S* X x(%/^p) + Sx{l- x(V7V'p)) , 



yields the result. 



□ 



^Here the estimates become (5-dependent, because a large 5 requires a large A'^. It is 
uniformly controlled when 5 < 5q. 
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2.3 Unitaries and effective Hamiltonian 

We strengthen a httle bit the assumptions (|2.4p - (|2.5p . with the condition 

no(g, T, e) = uoiq, T, e)P+uo{q, r, e)* (2.23) 



with 



fulfiUed in our example. The operator uq is nothing but the local unitary trans- 
formation Uoiq, t)^ on L'^iW^; C^) . 

With the approximate projection 11 ~ Il{q,eDq,T,e) given in Proposition 12. 3( 
Proposition IA.5I tells us that a true orthogonal projection P can be associated 
when So is chosen small enough, by taking 



1 
2hT 



\z-l\ = l/2 



(z-U)-^ dz, 



with 



P = P((7,£D„T,£) , Pe5„(l,g,;M2(C)), 
P(g, p, r, e) - n(g, p, r, e) = e^+'^Pi (g, p, r, e) + e^+^s^^ (^^ ^) ^ 

^ ,g,;A^2(C)), 



(2.24) 



(2.25) 
(2.26) 



with 



i?i , i?2 G Su 

Pop = P = P* , 



H,P 



-3+2(5 



with 
and 



CuC2 e Su{hgr;M2{C))., 



P - uoP+Uq 



< Ce. 



(2.27) 



(2.28) 



For a general fi G C^, R2 and C2 are included in the main remainder term. 
When /i is quadratic in {\p\ < r}, then one can assume that Ri and Ci vanishes 

m for r' < r, according to Proposition l2.4l and Proposition 

Instead of constructing unitaries between P(t, e) and P+ by the induction pre- 
sented in [PST] and similar to ([2J3)) . ((2111) . ((2J3)) . ((2J6)) . we use like in jMaSoj 
Nagy's formula f fN^ . fMi55] . [FST] ) 

P2 = WPiW* , W*W = WW* = 1 , 
with W^{l-{P2-Pif r^^^P2Pi + {l-P2){l-Pi)] , (2.29) 
when P,j^P'j = P* forj = l,2, and jjPj - Pij|£ < 1 , 

easier to handle for direct second order computations in our case. 
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Proposition 2.5. With the definitions (|2.23p and (|2.24p a/ter Proposition 
there exists a unitary operator U on L^(R;C^) such that 

P = UP+U* , 

U ^U{q,eDq,T,e) , e 5„(l,.9r;A^2(C)) , 
wi/i U{q,p,T,e) = ua{q,T,e) + sui{q,p,T,£) + e^U2((7,p, r, e) , 

uiiq,p,e) = f, (5^^:no)Mo, 
and e"^*\ii,e~^'^U2 G S'„[ — = — ^ ,gr;A^2(C)] . 

{J^){V^p)°" 

Moreover, when fi is a quadratic function in {\p\ < r} and r' is fixed in (0,r), 
the term U2 can he decomposed into 

U2 {q, p,T,e) = e^^v2 {q, p, T, e) + e^^v2 {q,p,T,e) 

where V2 does not depend on p in and V2 , V2 belong to the 

symbol class Su ( — — ^ , ; 2 (C) ) . 

Proof: The notation R will denote a generic remainder term of the form R = 



R(q, eD„, r, e) with i? e 5„ — j= — , g^-; A^2(C) . The notation R is 

used for a symbol R, like R but which vanishes around 

apply Nagy's formula (|2.29p with Pi = uqP+Uq = Uo{q,T,e) and P2 = P{T,e) 
with 

P = Uoiq, T, e) + elli {q, eDq, r, e) + e^Ha {q, eD^, r, e) + e^+^^R + e^+^^R . 
In the expression of VF(e) given by (|2.29p . the first factor is nothing but 

(1 - (P - no(g, r, = 1 + ^(ni)2(q, r, e) + e^+'' R , 

owing to P - Ho = Osle^*) • In the factor [Pa^'i + (1 - -P2)(l - ^"1)], the first 
term equals 

P o Ho {q, e) = Ho {q, e) + eHi {q, t, eZ?,, r, e) o Hq (q, r, e) 

+ e^Hs (g, eD„ r, e) o U^iq, r, e) + e^+^^P + e^^+^^P , 

while the second term is 

(1 - P) o (1 - no(<z, T, e)) = (1 - no(g, T, £)) - elii (q, eD^^T, e) o (1 - no(g, r, e)) 
- £2n2(q, ei^^, r, £) o (1 - no((Z, re)) + e^'+^^R + e^+^^R . 
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Hence we get 

W ^ [P2P1 + {I - P2){1 - Pi)] = I + eIli{q,eDg,T,e) o asiq^T) 

+ e^n2{q, eDg, r, e) o asiq, r) + e^+^^R + e^+'^^R . 

The operator U{e) is given by U{e) = W{e) o uq . The semiclassical calculus 
recalled in ()A.2p - (|A.3|) yields the result. In the decomposition of U2, when /i 
is quadratic around 0, the terms which are linear in e^* come with the second 
derivative of /g, which does not depend on p . □ 



Proposition 2.6. Introduce the notation for k G {!,...,£?} 



Ah Xh 



IfU is the unitary operator introduced in Provosition \2.5\ the conjugated Hamil- 
tonian U {e)* H {e)U (e) equals 



U*HU = 



h+ 
/i_ 



e'+^'Ri + e^+^'R2 



(2.30) 



where the remainder terms Ri_2 = Ri.2{q, ^Dq, t, e) belong to OpSu (1, gr', A^2(C)) 
and additionally Ri = in when fi is quadratic in {\p\ < r}, 

with r' < r . 

The symbol /i+ and /i_ are given by 



Up) + E+{q) - e{dj,fe).A + -{d'^^^J,)AkA, 
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^'(dlpje)^^ , e'id,,.fe)idM 



-XkXi 



E+ -E. 



XkXi, (2.31) 



h^ = fM + E^{q) + e{dph).A + {dl^^J,)AkA, 



-AfeAf AfcAf . yi.6l) 



Proof: From the semiclassical calculus, we already know that U*HU is 
a semiclassical operator with a symbol in Su{{^J^q') , gr', ■M2{C)) . Its off- 
diagonal part equals 



[l- P+)U*HUP+ + P+U*HU{1- P+) 



U 



P, 



P.H 



U. 



The almost diagonal form (|2.30p of U*HU is then a consequence of (|2.27p . 
For the second result, it is necessary to compute the diagonal part of the symbol 
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U*fHfU up to C'(e^+^'^) in Su{l, gr] M2) ■ Let us compute the diagonal part 
of 

B := {u* + eul + e^u*^)fHf{ua + eui + e^U2) , 
or cquivalcntly (uqBuq)^ with our notations. 

Since Uq = uo{q,T,e) and H = e^^ J\{\J t' t" -p) + V((7, r, e) , the first Moyal 
product equals according to (jA.2P - (jA.3p . 



.3+25 



R 



UqH + £ 



1 



2i 



where i? and R denote generic element of Su{^^gr\ A^2(C)), with the additional 

property that R vanishes in when /i is quadratic in {\p\ < r} 

with r' < r. The reason for the possible decomposition of the remainder, comes 
again from the fact that the third order remainder term, proportional to e^'^'^^ , 
arises with the third derivative of /g, the second derivative w.r.t p of ui and the 
first derivative w.r.t p of it2 . 

In the same way, the complete expression of B is given by 



B{e) = UqHuo + e 



2i 



-u\Hui - ^-^^^{dqkuDui + Y- {u*i,H}uo 



By recalling that (uiUq)^ = by Proposition 12. 5[ we get 



,3+25^^^3+45^^ 



D 



z^B?+e^+'''R 



-3+45 



R. 



where B2 is made of several terms to be analyzed. We need the relations 

dqjul = -ul{dqjUo)ul , (2.33) 
^Ijqi'K = Kidq,'Uo)ul{dqmo)ul + ul{dqj uo)ul{d u^)ul - u^d^^^ Uo)u* 
= -idq,'K)i9q^uo)ul - {dqjU*){dy,'Uo)ul - u^d^^ .uo)u* , (2.34) 



and 



[uo{dqU*)] 



OD 



^{dqTlo)(T3 : 



(2.35) 
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coming from UqUq = 1 and the differentiation of UqUqUq ~ . 
For example, the first one simphfies the C'(£)-term into 

Many cancellations appear after assembling all the terms in B2 ■ We need 
accurate expressions for all of them: 

• By using again (uiUq)^ = 0, the term [UQU2H + Hu2Uq + uouIHuiUq]^ 
equals 

if, + E+)no{uou*2 + U2U*)no + if, + E-){1- Ilo){uou*2 + U2U*)il - Ho) 
+ (A + E^)Uo{uouluiu*)no + (/e + E+){1 - Uo){uouluiu*){l - Hq) . 

In the relation 

UqU2 + U2U0 + ului + ^ {uq,ui} + ^ {ul,uo} = £^+^'^i?i + e^+'^^Ri , 

At At 



the remainder terms satisfy Ri,Ri £ Su[ — 1^ — ,(7t-;A^2(C) 

with the same convention as for _R, R . This identity is obtained by writing 
that the 0{e^) remainder of — 1 vanishes and by noticing that the 

remainder i?i + e^^Ri involves second derivatives of ui w.r.t p and first 
derivatives of U2 w.r.t p . We obtain 

[U0U2H + Hu2Uq + uquIHuiu'^]^ = {E^ - E+){uquIuiUq)^ as 

- [uo {u*,ui}u* + uo {ul,uo} ulf + e^+^^R + e^+^^R . 



with R,R£ A^2(C)), again with the same convention. 

Again w 
factor is 



Again with {dpUi)uQ — {dpUiu'^)'-"^ ~ i(E^^-E ) '^q'' ^0 and p.35p . the last 



(52 j-^-j 



i{E+ - E- 

= -2^^^^^^(a,.no)(9,.no)a3 . (2.36) 
With uiUq = .^(g^'L^^ ) {dqkllo) , we have proved 

[uoU*2H + Hu2ul + UoulHu^ulf = ,no)(9,.no)a3 
- i/^^^(9,.no)(9,.no)a3 + s'+^'R, + s^+^' R, . (2.37) 
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• The term 



— UO {'«o/e''"l}"o + ^Wo{Ui/e,Uo}Uo "o(Og'= "oi^^l^O 

equals 

^Je [uo{Uq,Ui}Uq + Uo{uI,Uq}Uq]^ 

The diagonal part of (9^fcUi)wg is 

But differentiating the relation (9^fcno)no + Ilo{dqiIlo) = dgkHo w.r.t 
leads to 

(5'.,.no)^ = -{dg.nodgiiio + a,.noa,.no)(T3 . (2.38) 

\OD 9p^fe 



With (12351) and [uiuD ^ = i^^fz^^^f Hq, wc obtain 

[(a,. Ho) (a,. Ho) + (a,.no)(a,^no)] aa . 



i{E+ - 
We have found 



(2.39) 

The diagonal part of [{uqV, ui} + {u^V, uq}] Uq equals 



By using (|2.35p with 



(5,V)«^ = {E+ - E^){d,Uo) , 
and i;_no + ^+(1 -Ho) = V + (£:_ - £;+)cr3: 
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it becomes 

(a,. Ho) (a,. Ho) + [uo {u*,ui} u* + uo {ul,uo} u*of 

+ ^{E^ - E+) [uq{uI,Ui}uI + Uq{u{,Uq}uI]° (73 , 

The relation ([OSI yields 



1 d fe 
— [uo {ulV,Ul}ul+Uo {ulV,Uo}ul] = P'P' ^ {dqkIlQ){dql 



+ v-|?^^^(a,.no)(9,.no). (2.40) 



The term ^ [uq {u\, H}]^ is the sum of two terms 

^.[uo{ul,f^}f + ^.[uQ{ul,V]f . 

Since uodp^ul — (^e+-e^ off-diagonal while 



the second term equals 



The first term a priori contains more terms because Uq has to be differen- 
tiated: 

The first part is off-diagonal and vanishes after taking the diagonal part. 
By using again p.35p and ()2.38p and dqllo = {dqllo)'^^ , we obtain 

— -E, (CgfcllojlCqaiojcrs H (CgfciiojcTslCg^ffo) = 0. 

— h/— nj^ — hi— 

Hence we have proved 

l[u,{uiH)f = ^^(a,.no)(a,mo). (2.4i) 
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• The last term is 

T=i 2uo{idpJe)idqkUQ),uo}uQ 

Forgetting the i factor, the first part equals 2{dp^p^fs)uQ{dgkUQ){dqtUo)uQ . 
By using (|2.34p in the second part gives 

8T = 4:{dp^pje) [uo{dqkUQ){dqiUo)u'^]^ . 
Owing to (EH) and ([235)) . this gives 

With uoldgkUQ) ~ —{dqkU{))uQ, the last term equals 

(2.42) 

By summing ([07)) . (P^Bg)) . (PTiO)) . (piiTj) . (pii^ . we obtain 



92 f 

'^PkPeJ±, 



2^ (a,. Bo) (9,. Bo) + s'+^'Ri + e^+^'h ■ 



Bence the diagonal symbol that we seek, is 

'h+ 0\ ^ ff,ip,T)+E+{q,T) 
h.J { f,{p,T)+E.{q,T)^ 

D fe D D 

- ie{dp^fe)ul[{dqkUa)uo] ito - e^^^^|^Uo [(9qfcUo)itS] [{9qkUo)u^] uq 

+ £'%^^^%^«5(9,.Bo)(a,.Bo)a3Ko + s'^^u^id,kUo){d,mo)uo . 

The symbol [(5qfcUo)uo] ^ equals 

wSno(9^fcWo)it5BoUo + Uo(f - no)(9qfcUo)w5(f - Bo)wo 

= p^ui{dqku^)p^ + (f - p+)u*(9,.wo)(i - ^+) = (^0' _y • 

For the last term we deduce from Bo = uqP+Uq and {dqUo)uQ + uo{dqUo) = 0, 
P+ [u*Q{dqkIlo)uo] [ul{dqtIla)uo] P+ 

= P+ [K{9q>'^Q)P+ + P+{dq''K)^Q\ [K{9q'^a)P+ + P+{dqeul)uQ] P+ 
= P+ul{dqkUt^)P+UQ{dqiUo)P+ + P+U*Q{dqkUo)P+{dqlul)uoP+ 
+P+{dqkul)U(iUl{dqtU(i)P+ + P+{dqkU*^)UaP+{dqlul)uQP+ 

= -P+ul{dqkUo){l - P+)ul{dqtuo)P+ = XkXg . 
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Taking the bracket with (1 — P_|_) is even simpler and gives 

(1 - p+) [uUd,^no)uo] [u*idgmo)uo] (1 - p+) 

= -(1 - P+)u*{dg.uo)P+u*{dgeuo){l - P+) = Me . 
This ends the proof. □ 



2.4 Discussion about the adiabatic approximation of the 
Born-Oppenheimer Hamiltonian 

The Theorem 12 .11 is a direct application of Proposition 12.61 by taking 



feip) 



_2(5 



r'T"\p\Hr'r"\p\' 



The operators 



h±,Bo{q,£Dq,e) 



-2<5„ 



.k=l 



Weyl 



fe=i 



(-v).A 



A.(-V) 



\x\ 



is nothing but the usual adiabatic effective Hamiltonians which can be found in 
the physics literature, including the Born- Huang potential \X\'^ ~ Ylt=i l^feP ■ 
We refer to [PST| and |PST2j for a discussion of the various presentations of 
the calculations and additional references. 

Even in the region |Vt't"|p| < r^j with p quantized into eDq, this approxi- 
mation makes sense, only for 5 > Q, because of the additional term 



T2 



XkXt{edq^){edq,) 



coming from the last terms of ()2.3ip and p.32p . that we have included in the 
remainder. It is not surprising (see |Sorj . jMaSoj ) that the degree of the dif- 
ferential operators increases with the degree in e in the adiabatic expansion of 
Schrodingcr type Hamiltonians. The argument of physicists says that this ef- 
fective Hamiltonian is used for relatively small frequencies (or momentum) so 
that XkXip'^p^ is negligible w.r.t \p — eA\^ + e^jXp . The introduction of the 
additional factor e^^ with 5 > provides a mathematically accurate and rather 
flexible implementation of this approximation. 



3 Adaptation of the adiabatic asymptotics to 
the full nonlinear minimization problem 

In this section, we adapt our rather general adiabatic result to our nonlinear 
problem. In a first step, we give an explicit form of Theorem 12.11 in our spe- 
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cific framework. Those results are effective when applied with wave functions 
localized in the frequency variable, ip = x{y/'''x'''ySDq)ip for some compactly 
supported x ■ It could suffice if we considered minimizing the energy among 
such well prepared quantum states. We can do better by using a partition of 
unity in the frequency variable, which will be combined, in the end, with the a 
priori estimates coming from the complete and reduced minimization problems. 
Finally an estimate of the effect of the unitary transform U on the nonlinear 
term is provided. 

3.1 Adiabatic approximation for the explicit Schrodinger 
Hamiltonian 

Let us specify the result of Theorem 12.11 by going back to the coordinates 
q = ((?', q") — {x, y) and r = (r', r") = [r^.Ty) . Provided that 14, r (a;, y) fulfills 
the proper assumptions, the unitary transform introduced in Theorem I 2 . 1 1 trans- 
forms given by (|1.9p into the Born-Oppenheimer Hamiltonian with a good 
accuracy in the low frequency region, that is when applied to wave functions '0 
such that = x{y/V^£Dq)ip for some compactly supported x ■ 
The operator is the e-quantization of the symbol 

-25, , i„i2 , „, / ^ ( E+{q,T,e) 



e TxTy\p\ +UQ{q,T)\ Q ' £;_(g^r,p^ ' ' 



with E±{q,T,e) Ve^rix.y) ± n{. —x) = Ve^r{x,y) ± . 1 + — . 



The operator 7/0(9, r) equals 

Uo{q,T) = (^g^i^ ^^fj ) =Uo{q,T)* 

t 
with C = cos(-) , S = sin(-) 




Proposition 3.1. Assume 5 G (0, (5o], {Tx,Ty) E (0,1]^ and assume that V^^r 
belongs to S'„((. /^x), — - + ^dy^) then the matricAal potential 

,w ^ . \ , nf \ f cos{9) e"^sin(6')\ 

V(g,r,e) = V,Ax,y) + n{^-x) (^^-., _ ^^^^^^^ ) 



fulfills the assumption of Theorem \2.1\ 

Choose the function 7 G (R) such that j = 1 in a neighborhood of [— , r^] , as 
in Theorem ] 2. 1\ and consider a cut-off function x G C5"((— r^, r^)). When U = 
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U{q,eDq,T,e) G Op5'u(l, f/r;-A^2(C)) is given in Theorem \2.1\ the identities 



.2+25 







Zjl y. ^ 1 IJL y. 



\x{TxTy\eD,\^) , (3.1) 



hold with 







e V 2 u_ g \ 



:^ \ U*x{TxTy\eD,\^),{i.2) 



-dl -\dy+ l- 



2jl 



'x'y 



-dl - a, 



K,r(x,y)-Jl + ^X2 



and the estimates 



r^(l + ^x2)2 r,(l + l-a;2) ' 



l|i?i(r,e)|U(L2) + ||i?2(r,e)||£(i2)<C 
which are uniform w.r.t 5 G (0, ^o]; G (0, 1]^ and e G (0, Eq] • 

Proof: Following the approach of Section [2J the Hamiltonian Hun is decom- 
posed into 



H 



Lin 



H{q, eDg, T, e) + R^{eDg, t, e) 



with H{q,p,e) = £ T^Tj^bl 7(7-xTy|pr) + 



uo 



E+ 




and R^{p,T,e) = e^^T^Ty\p\'''{l ~ j{T^Ty\p\''')) 



The Hamiltonian H{q,eDg,T,e) fulfills the assumptions of Theorem 12.11 with 



the metric g-r 



because we assumed 



T,, 



Is. t\ 2 T-r, 



^dy') 
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while the gap E+{q,T,e) — E^{q,T,e) equals 2y^l + ^x'^ . Actually the esti- 
mates 



-|a| 



are due to d^-O = — , I'f^T^ and 9„e"'' — i./^e"^ . Moreover the explicit 

computation with a -\- j3 — 1 leads to 



iul{dxU{)) 



idx 







and 



Ay^ Xy 

Xy -A 



) = iuQ{dyUo) = \ —S ( 
yj V V~ 



5 -Ce"^ 
-Ce-''^ -S 



(3.3) 
(3.4) 



Hence the effective Hamiltonians hBo,±iQTsDq,e), when restricted to the region 
{^Ta;Tj,|p| < r^}, are given by the symbols 



With cos(6i) 



pI + iPy T ej^ sin^(^))^ + i(|a.^p + ^ sin'm 



and sin(0) 



the Schrodinger Hamilto- 



nian corresponding to the RHS is e'^'^ TxTyC *27f??'^_|_g+*2;^y with 



H+ = -di- d„±i- 



2jl 



^ ye,r(a;,y)±Jl + - 



The remainder term in Theorem 12.11 is 

£2+4*i?i(g, eDq, T, e) + e3+2^i?2(<7, eDq, r, e) , 

with i?i^2 G S'„(1,5t; ■A/i2(C)) and where i?2 vanishes in a neighborhood of 
IpI < r^} . The first term provides the expected 0{£'^^^^) estimate in 
L^(]R2;C2). 

It remains to check the effect of truncations. All the factors, including the left 
terms 



pI + (pj,Tesin2(-))2 



(l-7(r.r,|pp)) 



belong to OpSu{{y/TxTypy , gr', M2{Cj) . For any a,b belonging to the class 
OpSu{{y/TxTyp)'^ , g-r; A^2(C)), where b vanishes in a neighborhood of | y/TxTy \p\ < r^}, 
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and two cut-ofF functions Xii X2 S ^^(("'^71 ''7)) such that xi ^ X2 (see Defini- 
tion |X]3|, the pseudo-differential calculus says 

(1 - X2(r.rjpn)f af Xi(T,T^|pn G Nu^g^ , 
bfxi{rxTy\p\'')&Uu^g^ , 

with uniform estimates of all the seminorms w.r.t r € (0,1]^ and S € (0,So] . 
Applying this with Xi = X and various X2 such that Xi ^ X2 ^ 7i implies that 
the remainder terms due to truncations are 0{e^) elements of C{L^) for any 
N gN, uniformly w.r.t r G (0, 1]^ and S € (0, 60] . Fixing N>2 + 4So ends the 
proof of (IXTj) . 

For p.2p use p.ip with a cut-off function xi such that X ^ Xi and conjugate 
with [/ : 



HL^nUXl{TxTy\eDg\^)U* 









Xi(T.r,|£i?,|2)J7* 



Right-composing with xiTxTyleDql"^) and noticing that 



XiiTxTy\sD,\'')U*x{rxTy\eD,\'') - C/*x(r:.Tjei?,n = ei?,, r, e) , 
with R E Mu.g^ lead to p.2p like above. 



□ 



3.2 Linear energy estimates for non truncated states 
Proposition 3.2. Assume 6 G iO,5o], t ~ {Tx,Ty) e (0,1]^ and that V^^r be- 



longs to the parametric symbol class Su{{a/^x), 



^dy^). Setx 



xiTxTy\eDq\'^) for x G r^, r^)) . When U is the unitary semiclassi- 

cal operator U{q,eDq,T,e) £ OpS'„(l, A^2(C)), given in Theorem \2.1\ and 
parametrized by a truncation in {^TxTy\p\ < r^} , then for any x € C^((— r^, r^)) 
the estimates 



^, [U*Hi^,,,U-s'+''rxryHBo]4') < Ce'+'' [e'+''\mh 

+ 11(1 - m\h +!l(l - m\L^ X \\^\eDq\{l - xM 
V^, [Hun-e^+'^^TxTyUHsoU*]^ 



< Ce 



1+25 rgl+25 



\h 



(3.5) 



(3.6) 



+ 11(1 -X)^lli2 + 11(1 -X)'^IU^ X \\^\sDq\il-x)^\\j 
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hold uniformly w.r.t r E (0, 1]^ and S G (0, Sq], for all ip G L^(M^; C^), with 



H 



BO — 



e"'^^i?_e*^^y 



H± = -di- \ dy±i 



2 



Proof: Set 



V e"'^^^i7_e'^^^. 



and bound the terms {xip , "Dip) and ((1 — xip , Vxip) by Ce2+4'' ||-0||2 ^ with the 
help of Proposition 13.11 The remaining term is 

{{l-x)i>,V{l-x)i,), with x = x{r.Ty\eDq\''). 

The operator V can be decomposed according to 2? = e^^^^TxTyVkm + 'Dpot 
with 

^^'"^ - ^ l^d |i?, + A|2 + |X|2j ' ^•^■^^ 



2?pot - U*Uo U*U ~ VeA^,y) ^1 + ^x2 (^^ ■ (3-8) 

The normahzation in p.Sp allows to use directly the semiclassical calculus if 
one remembers that UqU = Id+sR and (j*Uo = Id+e^' with £~'^^ R, e^'^^ R' e 
S„(-7^,5r;X2(C)) and E±{q, r,£) = ± ./l + ^^2 . We obtain 

i((i - x)^.v,ot{i - x)m < c£^+''ii(i - m\h , 

which corresponds to the second term of our right-hand sides. 
The kinetic energy term (|3.7p is decomposed into 

with 



'Din = e^'U*T,Ty\eD,\^U~e^'u;T,Ty\eD,\^Uo, (3.9) 

'\D,-A\^ + \X\^ 



(3.10) 
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Writing p.9p in the form 

T^ln - e^Hu* - US)T,Ty\eD,\^U + e^'UST,Ty\eD,\'{U - Uo) , 

while T,Ty\eDg\^ e Op5„((7^p)2,g,;Al2(C)), [/, J7o e Op5„(l,.9,;M2(C)) 
and £-i-2i(f)_ ^ OpSJ , .,\ ^^ ,gr;M2{C)) , leads to 



|((i - x)^, - x)V^)| < - x)V'l!i. , 

which is even smaller than the 'Dpot upper bound. 
In p.lOp . the first term can be computed via 

($, U^\D,\^Uo{qm = {Dg,{uoiq,T)<S>),Dg^iuo{q,T)^f)) , 

with Dg.{uof) = UQ{Dg^ - iu^dq-uo)/ and equals 

9je{2;,a} 

+ {iu*Q{q)dq^u„{q)y 
Meanwhile expanding the entries of the second term in p.lOp gives 

\Dq T A{q)f = [^l T A, {q)Dq^ T Dq^ A, {q) + A, {qf . 

By using the expressions p.3p and p.4p for and iugOqUo, we obtain 

/cm 



2 Vi?+ 

with i?± = ±i(i:)^(a:,6l) + {dJ)D.^)e^''P - sm{9){Dye^''^ + e^'^Dy) 



and (c?^;^) = /• ' ^ — 5- , sin(6') 



Therefore, we obtain 

|((1 - x)V', - X)^)\ < 4max( /^)|l|D,|(l-;e)V'IUHI(l-x)V'l|L^ 

and, owing to Tr, G (0, 1], 

|((l-xK.,£2+2*r,r,I?L(l-x)^)| 



< 
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This ends the proof of 
For p.6|) it suffices to replace tp in (|3.5p by U*ip with a kinetic energy cut-off 
function xi such that x ^ Xi ^'iid then to use (1 — Xi)U*x G OpAfu.g^, with 
uniform scminorm estimates w.r.t r G (0, 1]^ and S € (0, Jq] • The L^-norm 
of the corresponding additional error term is 0{e^), for any A^. and one fixes 
TV > 2 + 4(5o . □ 



3.3 Control of the nonlinear term 

In this subsection, we estimate the effect of the operator U = U{q,eDq,T,e) 
belonging to OpSu{l, Qt', ■M2{'C)) on the nonlinear term dxdy . 



Proposition 3.3. Let U be the unitary operator introduced in Theorem \2.1\ 
The inequalities 

J mx,yt dxdy > {I - Ce'+^') J \{Ui^){x,y)\^ dxdy, (3.11) 
and J\il;{x,y)\'^dxdy>{l-Ce^+^^)J\{U*i}){x,y)\^dxdy (3.12) 
hold for any G L^{R^;C^) . 

Proof: For ipi and ip2 belonging to L^{M.'^;C'^), the local relations 

li'iWq) = {\Mq)\' + '2MM<i)Ai^i~i^2){q)) + \Mq)\'f 

= + 2|V^2|'|V'i ~ ^2|' +4 (^Me(^2 , ^1 - ^2) + IH'1 - V'2|') 

+41^-21' Me(V2, (V'l -^2)) 

> \MQ)\'~MMq)\'\Mq)-Mq)\, 

is integrated w.r.t q — {x,y) e M^, with Holder inequality, into 

llV^lllia- j dxdy>\\i'2tL.-^^2\\lA\i'l- i'2\\L^■ 

Wli'h ^1 = Uoip = uaiq)ip, \^iq)\^ \M<l)\^ for ah q G M^, and V2 = Ui' = 
tpi + {U — C/o)V', we obtain 

Mh = iiv^iiii^ > mwi^ - mi'iihuu-uomL^ . 

The operator U — Uq equals e^^^^r{q,eDq,T,e) with r belonging to the class 
Su( , ^ ,' ——,9r;M2{C)), where we recaU = + ^dy^ + 

(^'^J^^^pyx ■ After introducing the isometric transform on L^{^] C?) 

{Te,r'p){x, y) = e^TxTyip{e^T:rTyX, e^T^Tyy) , 
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the difference U — Uq becomes 



U-Uo = e'+%-^ri{eT,x,eTyy,D,,Dy,e,T)T,^, 

with ri uniformly bounded in 5(1, ^^r^ + d{Tyy)'^ + j^; 7W2(C)) . A fortiori, 

the symbol ri{eTxX, ETy-y, ^, ry; e, r) is uniformly bounded in ^(l, dq^ + ^pi) and 
the Lemma 13.41 below provides the uniform bound 



\\U-lfo\\ciL^)<Coe'+' ■ 

We have proved 

MU>\\mt.-ce'^''\\mumL*- 

which implies p.lip . The second inequality (|3.12l) is proved similarly with 
u* = U* + {U* - U*) . □ 

The result below is a particular case of the general bound, 1 < p < oo, for 
pseudodifferential operator in OpS{l,dq^ + C{'Hi;'H2)), Hi Hilbert spaces, 
stated in |Tay| -Proposition 5.7 and relying on Calderon-Zygmund analysis of 
singular integral operators. 

Lemma 3.4. For any p G (l,+oo), there exists a seminorm n on S{l,dq^ + 
j^;M2iC)) such that 

dn^ 

Va e S{l,dq^ + -^;M2(C)), \Hq, D,)\\c(l.) < n(a) . 

2 2 

Proof: The Proposition 5.7 of [Tayl says that for any a € S{l,dq'^+^ ;7W2(C)), 
the operator a{q, Dg) is bounded on L^'(M^; C^) . It is not difficult to follow the 
control of the constants in the previous pages of |Tay| in order to check that 
||a(g, Dq)\\c(Lp) is estimated by a seminorm of a . More efficiently, a linear map- 
ping from a Frcchet space into a Banach space is continuous as soon as it is 
bounded on bounded sets. Apply this argument with the result of |Tay[ to 

5(1, dq' + -^;M2{C)) 3a^ a{q, Dg) £ CiLP {R^ ; C^)) . 
KP) 

□ 



4 Reduced minimization problems 

In this section, we assume that the potential V^^r satisfies (|1.13p - (|1.14p . After 
the first paragraph of this section and in the rest of the paper, we focus on the 
case Ty = I, Tx Q (and e — > 0) . Two reduced problems have to be considered: 
1) the one obtained as £ — > and Tx is fixed; 2) the one derived from the previous 
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one as Tx ^ and which is parametrized only by (G, The hnear part 

of this latter reduced problem is a purely quadratic Schrodinger Hamiltonian 
(with a constant magnetic field), from which many a priori information can 
be obtained. This section is divided into three parts. First we specify the 
potential V^^r and check our main assumptions for the general theory. Then we 
review some properties of the reduced Gross-Pitaevskii problem parametrized 
by {G,iv) ■ Finally we make the comparison with the reduced Gross-Pitaevskii 
problem parametrized by (G, (ytTx) as t^: — >■ and deduce properties which will 
be necessary for the study of the complete minimization problem. 



4.1 Reduced minimization problems 

Lemma 4.1. The potential V^.r defined by (|1.13p - (|1.14[) belongs to the class 
Su (^{^J^x),gq^r^ with the metric gq^r = ^^([^"'1=;^^) + -^dy^ ■ 

Proof: After the change of variable (x', j/') — (y^x, ^J^y), it is equivalent 
to check 

_2 -^2 



.2+2(5 



-v{TyX, rj,J/) + Vl + x^-i 



.2+2(5 



(l + a;2)2 l + .x2 



eSU{x),^+dy^). 



It is done if v{TyX,T,y) G ^ + dy^) . Wc know v S 5(1, j^^) ■ 

Hence for all (a, /?) € there exists Ca^p > such that 



Vre (0,l]2,Va;,2/eR2, \d^d^ {v{TyX,T^y))\ < G„,^ 

1 



' X ' y 



(^+x2)"/2 



(1 + T^x^ + r2y2 

<G„,^(a;)i-l"l, 



which is what we seek. 



, + 13 



□ 



If the error terms of Proposition 13.21 and Proposition 13.31 are assumed to be 
negligible, the energy of a state ip ~ U* ( ^ ] is close to 



e'+'*r,r„(a_,ff_a_) 



Ge.T 



G 



J |a_|4 dxdy = e^+^^T,Ty£ria-) , 



with £r{a-) — (a_ ,H-a-) + ^ / 1^-1^ dxdy . 



(4.1) 



2J1 + 5=^x2 



+ v{^T^TyX, y/T^Tyy) . 

'-yTxTy 



with the potential v chosen from (|1.14p . 

When ^ and r^Ty are small, in particular in the regime 1 and = 1 that 
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we shall consider, this energy is well approximated by 



2 , x^ + y^ 



2 



+ (4.2) 



as this will be checked and specified in the next paragraph. Although more 
general asymptotics could be considered, we concentrate from now on the regime 

^ 1, Ty = 1 . The parameters Iv and G are assumed to be fixed as r^; -> . 

In order to prove that the ground states of £h and £e are close, we need 
good estimates on the energy £h- 

4.2 Properties of the harmonic approximation 

The energy functional £h does not any more depend on and is parametrized 
only by (G,^y) . Let us start with its properties. We introduce the spaces T-Li 
and 7^2 which are given by 

■Hs = \u^L\M') J2 lk"^9^"llL2 < +00 L .s = l,2,(g-(a-,y)) 

[ |q| + |/3|<s J 

(4.3) 

endowed with the norm = X]|q|+|^|<s Ik^^g^llia • For a compact set 

K of Hs and for u £ Hs, the distance ds{u,K) follows the usual definition 
min„gif ||m — . The self-ajoint operator associated with the linear part of 
£h is denoted by 

H,„..=-dl-{dy-'fr + ^^±^, (0<^v<+oo). 

Its domain is 7^2 while its form domain is Hi . Note also the compact em- 
beddings 7^2 CC Hi CC L-^ n L'* . Following the general scheme presented in 
|HiPrl |Sjo] , its spectrum equals 

aiHe,.) = {(1 + 2?7,+ )r+ + (1 + 2n_)r_ , K,n_) € N'} 



Proposition 4.2. The functional £h admits minima on {u e 'Hi, ||wj|i2 = 1} , 
with a minimum value £H.m,in satisfying 

^H.niin > r+ +r_ > ^ . 

The set of minimizers Argmin £h is a hounded subset of 'H2 and therefore a 
compact subset of {u € "Hi, ||u||^2 = 1} . Moreover for any ip € Argmin £h, f is 
an eigenvector of Hq + G|(^p . Finally there exist two constants C = Ci^,g > 
and V = I'lv.G G (0,1/2] such that the conditions u G T-Li, ||u||^2 = 1 and 

£h{u) < £H,rrnn + 1, imply 

dn, {u, Argmin £„) < C{£h {u) - f h,„„;„)" . (4.4) 
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Proof: On Hi and H2, the scalar products 

(u, v)i^e^, = {u, He^,v)L2 

provide norms k = 1,2, respectively equivalent to . In this 

proof, all the "uniform" estimates are actually parametrized by {G,lv)- The 
nonlinearity ^ / dx as well as the constraint ||w||£,2 = 1 are continuous 

functions on n L'^ while the quadratic part of Eh{u) is simply with 

£h{u) > \\u\\l^^^ > r+ + r_ > ^ . The compact embedding Hi CC n 

thus implies that the infimum inf„g-^j ^ ||u||^2=i £h{u) is achieved. 

A minimizer ip €! Argmin £h solves in a distributional sense the Euler-Lagrange 

equation 

where is the Lagrange multiplier associated with the constraint HiysHis — 1 . 
By taking the scalar product with Lp, one obtains the bounds for X^: 

^H^rain — X^ ^ H ,min • 

Since "Hi is also (compactly) embedded in L^(M-^), the equation 

Hi^.-P = -G\ip\'^ip + \ip 

ensures that ||iy9|j2.fv uniformly bounded on Argmin . Therefore Argmin Eh 
is a bounded subset of ^2 and a compact subset of . In an Hi -neighborhood 
of iy9 e Argmin En (||<y5j|L2 = 1), the L^-sphcrc {u G H\, ||u||i2 = 1} can be 
parametrized by 

Notice also that the potential G|</Jp is a relatively compact perturbation of , 
so that Hi^ + G|(^p is a self-adjoint operator in L^(IR^) with domain H2 and 
with a compact resolvent. With ((^s , v)\^i^ = — G Jjja I'/Sp^'^ dx for ip is an 
eigenvector of + GYp\^ and v _L p^ the energy Eh{u) becomes 

~2Gil~\\v\\%)Rc [ \p\^pvdx+^[ Kl-ll^lliO^ + H'rfa; 

where v lies in the closed subset Hi^^ ~ {v e Hi, {tp , v) 1^2 ~ 0} of "Hi and F^{v) 
is the composition of the compact embedding Hi — > L'^OL'^ with a real analytic, 
real-valued, functional on L'^ O . Hence on Hi,ip endowed with the scalar 
product ( , )l^e^,, the Hessian of Eh{{1 — \\v\\^l2)p + v) equals Id (0), 
with D^Fiy (0) compact (and self-adjoint) . We can apply the Lojasiewicz-Simon 
inequality which says that there exist two constants dp > £ (0, 1/2], such 
that 

\\v\\l4v < {£h{{1 - \\v\\l2)p + V)- EH.rmnY^ • 
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Since the set Argmin £h is a compact subset of T-Li, it can be covered by a finite 
number of neighborhoods of tf^ g Argmin £h, 1 < « < iV, where a Lojasiewicz- 
Simon inequality holds. Take 

vg., a = rnin v^. and Ci,, n — ^ max . 

l<i<N ^ ' l<i<N ^ 

□ 

Remark 4.3. The Lojasiewicz inequality is a classical result of real algebraic 
geometry (see a.e. \Loj\ \BCRf ) proved by Lojasiewicz after Tarski-Seidenberg 
Theorem. It is usually written as |V/(a;)| < C|/(a;)|'^ with v £ (0,1] for a real 
analytic function of x lying around xq with f{xo) = . The variational form is a 
variant of it. It was extended to the infinite dimensional case with applications to 
PDE's by L. Simon in \Simf . We refer the reader also to \Chi[ \HaJe[ \Huaf and 
]BDLM^ for recent texts and references concerned with the infinite dimensional 
case or the extension with o-minimal structures. 

The nonlinear Euler- Lagrange equation is usually studied after linearization via 
the Liapunov- Schmidt process. Here using some coordinate representation of 
the constraint submanifold, especially when it is a sphere for a simple norm, 
allows to use directly the standard result for the minimization of real analytic 
functionals. 

When the minimization problem is non degenerate at every ip G Argmin £h, 
i.e. in the present case when the kernel of Id +D^F^p{{)) is restricted to {0}, the 
compact set Argmin £h is made of a finite number of point. When £v is fixed 
so that r_|_ and r_ are rationally independent, the spectrum of H^^, is made of 
simple eigenvalues and when G is small enough, G < Giy , the non degeneracy 
assumption is satisfied via a perturbation argument from the case G = . For 
large G, we can only say that the set of {G,£v) S (0, +oo)^ such that all the 
minima are non degenerate, viv,G — 1/2, is a subanalytic subset of IS? . In 
our case with a linear part H^^, which is a complex operator with no rotational 
symmetry, no standard methods like in \AJE!I allow to reduce the minimization 
problem to some radial nonlinear ODE. 

From the information given by the Lowest-Landau-Level reduction, when G 
and ly are large, the supposed hexagonal symmetry, after removing some trivial 
rotational invariance, of the problem (see \ABN\ \Nief } suggests that there are 
presumably several minimizers. 

A change of variable (p{x, y)e~^^y/'^ = au{ax, ay) with = 1/{IvVG) leads 

to 

£Hi^) = -±=£H{u) = -^l \{V-'-£vVGe,x^-^u\^ + GirM' + ^\u\'), 

(4.5) 

with the notations r = ?- = \r\ . This implies that £vVG is equivalent to a 
rotation value. We have the following results from the literature 
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• when iy\fG is small and G is large, the minimizer is unique up to rotation 
and vortex-free [A JR| : namely u(x^y) = /(r)e"^ for some real number c, 
where / does not vanish. If iyVC = 0, this is an adaptation of a result 
of [BrOs] . When iyVC is non zero, this requires refined estimates for the 
jacobian. 

• when £yVG is large, then vortices are expected in the system and this 
can be analyzed in details in the LLL regime (lowest Landau level) if 
additionally VC/iy is smaU (SB lABNj . More precisely, if VC/iy is 
small, then 

inf - ^ - inf ^lll = o (^^^ (4.6) 

where 

Elll{u)^ J Gir^\u\^ + ^\uf) (4.7) 

for functions u such that u{x,y)e^^'^^ is a holomorphic function of 
X + iy. This space is called the LLL. If u is a ground state of £h and w its 
projection onto the LLL, then \u — w I tends to in and C°'" as VG/iy 
tends to 0. If additionally, iyVO is large, then one can estimate inf iJ^LL 



|ABN| thanks to test functions with vortices and inf E'lll — O 




• if iyVC is large, and VC/iy is large, then this is a Thomas Fermi regime 
where the energy can be estimated as well |Aftj and is of order \fG j ly . 

We complete the previous result with another comparison statement which 
will be useful in the sequel. 

Proposition 4.4. There exists C = Ci^^^g > such that when u € T-Li satisfy 
^niu) < £H,min + 1; = 1, and solves 

Hev,GU + G\ufu = XuU + r 

with A„ G M and r ^ , then 

• u GH2 ; 

• there exists uq G Argmin £h, with Lagrange multiplier Xuo, such that 

|A„ - A„J + ||m - Un\\-H2 < (||r||i2 + {£h{u) - SH^minY) , 

where v = i^e^.c € (0, ^] is the exponent given in Proposition \4.2[ 

Proof: Since Argmin £h is compact. Proposition 14.21 alreadv provides uq € 
Argmin £h such that 

\\u - Uo\\-Hi < C{£h{u) - £H,minY ■ 
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Taking the difference of the equation for u and the Euler-Lagrange equation for 
Mo, we obtain 

Hi^-,g{u - uo) = (A„ - Xu„)uo + Xu{u - uq) + G{\uo\'^uq - \u\'^u) + r. 
Taking the scalar product with uq, with 

|||uo|^ttO - < C {£h{uq) +£h{u)) \\u - Uo||«i < C {£h{u) - £H,m%nT 

imphes 

Using the cUipticity of H^^-^q and the equivalence of the norms ||<y3||-H2 ^-nd 
\\Ht^-,GV\\L'^ ends the proof. □ 

4.3 Comparison of the two reduced minimization prob- 
lems 

In the regime Ty = 1 and — > 0, while > and G > are fixed, we compare 
the two minimization problems for the energies £t- and £h defined in (|4.ip - (|4.2p . 
We start with the next Lemma which is a simple application of the so called 
IMS localization formula (sec a.e. |CFKS| ). We shall use the functional spaces 
defined by (|4.3p associated with £h as well as the standard Sobolev spaces 
H^{R^) associated with f^, with s = 1, 2 and C i/"(R2) . 

Lemma 4.5. Letxi,X2 G C|f'(R^) satisfy X1+X2 = 1; suppxi C {x'^ + < l} 
and take a £ (0, ^] . Then the following identity 




+ g[ {xlxlKr:.)\u\\ (4.8) 

holds for all u £ Hi, with the same formula for £t{u) when u G . 
Moreover, £r and £h satisfy 

£r{u)^£H{xi{r^-)u)+£r{x2{r^-)u)-Trj2 I \{^ Xj){r^ .)\M'' 

+ G [ ixlxl){r:.)\u\' + R{u), (4.9) 

for all u e iJi(R2) with 

\Riu)\ < \ (£r{xi{<.)uY'^+£H{xi{r:.)uf'^) h|U.ri-3a . 
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Proof: The first identity is a direct application of tlie IMS localization formula 
(see a.e. |CFKS| ) which conies from the identity 

Px^p - xP\ = [P. X? - \ [x\P] P-\p [p. x'] 

when P is a differential operator of order < 1 and x is a C°° function. Simply 
combine it with the identity 

\ut = \xi{r:.)u\' + \X2{r:.)u\' + 2x?X^(t,".)|«|^ . 

Using the same argument for £t provides the same identity after replacing £h 
with £r, and it suffices to compare fT(xi(T".)w) with £'/f(xi(T".)u) . The defi- 
nition (|1.14p of the potential v and the condition a < \ imply 

■"{tI/^.) = r,.(x-2 + y2) on suppxi(r^".) ■ 
Therefore, we obtain, by setting Ut ~ xi(r".)ii , 
\£r{xi{r^-)u)-£H{xi{<-)n)\ = 

\{dy-l 



<[\\{dy-l 



2^/T+T,rX■ 
X 



^)ut\\l^ + \\{dy~i-)ur\\ 



L2 



1 + Vl + TxX' 



rXi(r:.)"ll 



< 



\{£r{xi{r^.)uY'^+£H{xi{r^.)uf^ 



□ 

Proposition 4.6. For any given {iv,G) € (0, +03)^, there exists ti^,.g > 
such that the following properties hold when Tx < ti^^g ■ 

• The minimization problem 



admits a solution u G H (R ) 



inf £t{u) 

l('l^>2^ 



• A solution u £ iJ^(R^) to the above minimization problem, solves an 
Euler- Lagrange equation 



X 



I 1/2 ^ 
tyTx 



2^1 + TxX'' 

with < A„ < 2£r.min md belongs to 



G\u\ 



u = X„u 
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• Moreover the minimum value Er.min = niinue_ffi(R2) . ||u||^2=i ^•r(u) satis- 
fies the estimate 

• For u e Argmin Er and any pairs x ~ (Xi:X2) iri CJj"(M^)^ such that 
Xl+Xl 1 wi/isuppxi C +2/^ < 1} andxi = 1 m {x^ + y"^ < 1/2}, 
the functions Xji'i'x^^ ■)u, j — 1,2, satisfy 

\\X2ir'J'.)u\\h < C^,e.,Gr'J' (4.10) 

^r(Xl(^'/'-)w) < fr,mm + Q/v.GT^/^ (4.11) 

c?«2(Xi(Ti/^)u, Argmin £h) < C^u ^gtT''''' , (4.12) 
^^v,Ge(0,i]. (4.13) 

^ constant Ca.b,c is a constant which is fixed once (a, 6, c) are given. 

Proof: Fix ty and G. Wc drop the indices €y,G in the constants. The 

exponent a will be fixed to the value ^ within the proof. 

First step, upper bound for inf{fT-(it) , u e H^iM.^) , \\u\\l^ = 1}: 

Let X = (X1jX2) and x = (xijX2) be two pairs as in our statement such that 

Xi ^ Xi according to Definition IA.3I Take ug € Argmin £h d'Hi- According 

to Proposition 14. 2[ it belongs to a bounded set of 'H2 so that || Igl^Molli^, with 

q = {x, y), is uniformly bounded. Hence, ^ supp Vxj U suppxiX2 implies 

|Vx,(rr.)PKP = 0(rr) while / x?X^(r:.)KI^ > . 

Lemma [4751 above with the pair x and a G (0, \] gives: 

^H^rnin = £h{uo) > f H (xi (t" . )wo) + £h(X2 (t" .)uo) " Ct^"" . 

On suppx2(t".), the potential is bounded from below by ^775:7 . Thus 

we get 

1 



'H.,rmn (||Xl"o|li2 + |lX2Mo|li2) > f H,m™ || X^l^O jl la + 7;^^ ||X2"oll 12 - CtI 

and finally 

||X2^io||i2 < C"tI^ , ||xi"o||i2 = 1 + 0(r«") , 

as soon as < {C £H,rmn)~^^'^°' ■ 

The function ui ~ Hxi'^oll/^l Xi'"o is normalized with 

£h .min <£h{ui) <£H.rmn+0{T°^), 



6a 
X 1 
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and %i(r".)ui = ui, X2iT^.)ui = . Applying the second formula of Lemma l4.5l 
with, now, the pair x, leads to 

with R{U,) < + i£H,,mn + 0(T^))l/')r,l-3" . 

With the estimate ^ < ^Hanin < C, we deduce 

It's time to fix a to the value | so that r^" — t^"^" = and 

inf Zriu) < £r{ui) < £H,7mn + KT^^^ ■ 

Second step - Existence of a minimizer: Once the function ui G i/^(M^) has 
been constructed as above, consider < tq with £h min + i^Tr. < ^i— ■ The 
functional ^ 

is the sum of a convex strongly continuous functional (and therefore weakly 
continuous) on i/^(R^) and a negative functional 



z;(ry2.)-l u). 



Due to the compact support of v—1, it is also continuous w.r.t the weak topology 
on i/^(M^) . Out a minimizing sequence (u„)nGN*, extract a weakly converging 
subsequence in iJ^(R^) . The weak limit, Uqo, satisfies 

with II Woo II L2 < 1- The same convergence holds also for the energy £r{u) — 
2f}r^ ll^lli^i SO that actually ||woo||l2 = limfe^oo HwritllL^ = 1 and Uqo realizes 
the minimum of £t{u) under the constraint ||m||l2 = 1 . 

The Euler-Lagrange equation can thus be written, with the stated straightfor- 
ward consequences. 

Third step - a priori estimate for minimizers of £r- 

Let u e i/^(R^) satisfy ||u||l2 = 1 and frCu) = £T,min < £h + kt^^^ . Take 
two pairs x = (XI1X2) and x' = (xi:X2): like in our statement, and such that 
x'l ^ Xi- The identities (j4.9p for £r and (|4.9p in Lemma l475l provide 

£r..rmn>£AXlir'J'-)u)^CiT^J' 

£„>£H{xiir'J'-)u)~C^ri/^ 

1 



^2/3 
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The first line says 

E^iXlirl"" ■)u) < £r,,mn + C^T^^' < Sn^mn + ^tI'^ + C^tJ/^ < C'^ , 

which combined with the second line provides the uniform estimate 

Therefore xi(ri^^.)u is uniformly bounded in "Hi with respect to Tx . Consider 
now the Euler-Lagrange equation 



1/2 



2 x/1 + Txx' 



p2 



■G\u\ 



U — XuU 



and write its local version for u'l =■ the form 



< - {dy - '{r 



u[ = Xuu'i - G\ui\'^u[ 

+ fu + fu , 



(4.14) 



with 
and 



J U 

f 

J u 



'ix{ 



l)x:{rl'\)dyu. 



4 1 + TxX 



2 "1 



after setting ui = xi(tJ^^.)u. Both functions, ui and therefore m'j^ = 'x!i{'^x^'^ ■)ui 
are uniformly estimated in "Hi and therefore in H^{M?) . From the embedding 
H^{M?) C L^{M?), the term G\ui\'^u'^ is uniformly bounded in L'^{R^) . For the 

term the support condition suppxi(''"i^^-) C < Tx ^^^^ imply 
while the estimate 

\\fl\\<C^-Tl''<C';, 

is straightforward. Hence the right-hand side of (|4.14p is uniformly bounded in 
L^(R^) and we have proved 

\Wi{r'J'.)u\\u. < Cl, (4.15) 

for any good pair of cut-offs x' = (xi; X2) ■ 

Fourth step- accurate comparison of minimal energies: 

2 /3 

We already know fr.mm < £H,min + ktx and we want to check the reverse 
inequality. Consider a minimizer u of £r and take two pairs of cut-off x — 
(Xi,X2) and x' = (xi,X2), such that xi < x'l ■ The identity (gS]) for £r and 
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(|4.9p of Lemma 14.51 used with x imply 

fr.™„ > Erixiirl"" ■» + £,{x2{rl"' .)u) 



,=1 Jvi' 



(4.16) 



■fr.mm > (Xl ^ •)") 

-rl"jzl \{Vx,){rl".)\M' + R{u). (4.17) 

After setting u'-^ = xi('^^^^-)'"; get the bound 

|(Vx,)(ri/^)p|.p== / KZM^^^!l^||,|V,p<C,ri/^ 

1/9 



while we already know from the third step the bounds £»{xi{tx ■)u) < C- 
when • stands for r or H, which implies |-R(m)| < C^t^^^ ■ From (|4.16p . w 
deduce, as we did in the first step with the energy Eh, 

\\X2{rl".)u\\l^<C^rl'^ , \\Xi{tI'\)u\\^.^1 + 0{tI'^), 

while the second line implies 

^T.min _ 1 ^ X } ^H,min — ^H,mi7i ^x x 

Fifth step- accurate comparison of minimizers: 

The function ui = Wxii'^x^'^ ■)u\\'^lxi{'^x^^ ■)u, satisfies 

while xi(Ti^^-)" solves the equation (|4.14|) for some pair x = (XI7X2) such that 

placing (x', u'^) w: 
setting ui = ||x('''^^^-)"llL2^Xi('ri^^-)" it becomes 



Xi -< Xi after replacing (x',m']^) with (x? Xi('nr^^-)''^) ■ After normalization by 



H.^^GUx + Gluipui = A„ui + + !l + /3 , (4.18) 
u,^\\x{rll'')u\\llxi{rl^'')u , ^ii ||x(ry9.)«||^i;fi(ry9> , 

= -2rV^(Vxi)(ry^).Vi.i - rl/\l:.xi){rl" .)u, , 
and /„3 = G||xi(ry^)^||i.(|Sip - l^i^ui + G(l - ||xi(ry^)^||iO|^.i| V ■ 

The estimate (|4.15p . for any new good pair x' = (x'iiX2) such that xi ^ 
Xi ^ xii implies that the terms and of the right-hand side of (|4.18p have 
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an L^-norm of order r^?^^ . For the third term the estimate (|4.15p also imphes 
that 

- l^^ilV = (l - , 

has an L^-norm of order O0^) (use the L°° bound for with \\\q\^ 

C^). We conclude by applying Proposition 14.41 . □ 

We end this section with a comparison property similar to Proposition 14.41 

Proposition 4.7. Let £v,G be fixed positive numbers and take u E i/^(M^) 
such that 

and which solves 

I 1/2 N 

-din - {dy - ^—±==f + ^^^^-^u + G\u\'u = X.^u + r„ 

with A„ G M and ||7'„||^2 < 1 . Then for any pair x = (xi, X2) G Cj^'{M.'^) so that 
Xi + X2 — 1 suppxi C {x^ + < 1} and xi = 1 m {x^ + < 1/2} , 
there exists T^^e^,G o,nd C^,i^ such that 

\\X2ir'J'.)u\\h < C^^iv.GT'J' (4.19) 
|£,(xi(ry^)?/) - £H.mn\ < C^j^^gtT (4-20) 
rf«2(xi(ry^)H, Argmin £n) < CxA.G(rf + llrulU^) , (4.21) 

when Tx < T^,iv,G O'nd where Viv,G G (Oj 5] exponent given in Proposi- 

tion\ 



Proof: The analysis follows essentially the same line as the study of the 
minimizers of £t in the proof of Proposition 14.61 By taking two pairs x' = 
(X1JX2) ^iid X ~ (XI1X2) such that Xi ^XIt^^ obtain successively like in the 
Third Step in the proof of Proposition [ 



£,{xi{rl'^ ■)u) + £r{x2{rl"' < Q ; 

-dl - {dy - f )2 + ^] u[ = X,u[ - G\ui\^u[ + fl + fl + ru where 
u'i,ui, /^'^ have the same expressions as in (|4.14p : 

• IIXi('^a^^^-)^ll«2 < owing to ||r„j|L2 < 1. 

From the last estimate, the refined comparison of energies like in the Fourth Step 
gives for a pair x = (xi:X2) such that xi < x'l' 

\\X2{tI" .)u\\h < G^tII^ £hM < £H,rmn + G^T^J' , 

with ui = \\xi{tx^^ ■)u\\~^xi{'''x^^ ■)'U' and G^Tx^^ < 1 for Tx < • The equa- 
tion ()4.18p is replaced by 

Hi^,GUi + G|ui| V = Xuui + fl + fl + fl + ||xi(Ty^)«|rV„ 
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without changing the expressions of /^'^''^ . Again the estimate ||Xi('''^''^-)"ll«2 ^ 
C^i is used with various cut-offs x'^, in order to get ||/^ + /^ + /^||l2 = 0{tI/^) . 
We conclude with the help of Proposition 14.41 applied to iti . □ 



5 Analysis of the complete minimization prob- 
lem 

We consider the complete minimization problem for the energy 



Ge.T 



1,4 I _ ^2+2(5 



and compare its solutions to the minimization of the reduced energies £r and 
£h, introduced in the previous sections. We work with Tj, = 1, — > 0, e — > 0, 
while £y, G and 5 E (0, Sq] are fixed. The analysis follows the same lines as the 
proof of Proposition 14.61 



5.1 Upper bound for inf {Ss^^i/j), = 1} 

The potential is chosen according to (|1.13p - (|1.14p while £v and G are fixed. 
The parameter is assumed to be smaller than t^^^q so that the minimal energy 
£T,min{Tx) oi Sr is achicvcd (see Proposition g^l) and \£T,min{Tx) - £H.,min\ < 
Ciy^G^x^^ ■ Moreover Proposition 14.61 also says that by truncating an element 
of Argmin £t, one can find a_ G "^2 such that 

||a_|U2 = 1 , |f,(a_) - £H,rmn\ < Ci^^gtT and \\a.\\n^ < Ce^^a ■ (5.1) 

Proposition 5.1. Under the above assumptions, take ^ U ( ^ v ) 

where a_ satisfies (|5.ip and U ~ U{q^eDq,T,e) is the unitary operator intro- 
duced in Theorem \2.1\ The estimate 

\£,i^P) - e'+''rx£r{a^)\ < C,,,g£'+"' • (5.2) 
hold uniformly w.r.t G (0, t^^.g] and S G (0, Sq] . 
Proof: Let us compare first the linear part by estimating 

By Proposition 13. 2[ it suffices to estimate 

1|(1-x)^||l^ and \\^\sDg\il-x)iA\L- 



2^1 + TxX^ 



V ' X 
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for some given cut-ofF function x G C^(R) with % = xi^xS^D"^) ■ Notice 



Hence by using the functional calculus of y/eDq, we can say that there exists a 
cut-off x' S C^(— r^, r^), with x' = 1 around and x' ^ X such that 

e'^(l-x)'e"'W < (1-x')', 
and e'W(r,e2|i5,|2)(l - x)'e-'^ < 2(T,e2|Z5^|2)(i _ ^,)2 + 3(1 - x'f 

as soon as e < eo < 1, for a convenient choice of eg and . By using ||a|| ^2(^2-1 < 
Ciy^G, we deduce 



||(l-x)V'lli2 < \\{l-x')a-\\l^<C,,,Grle\ 



and 



\\{V^e\DS{l-x)i>\\l-^ < 2\\{y^e\D,\){l-x')a-\\h+2\\{l-x')a-\\h 
By Proposition 13.21 we obtain 



^2+45 _|_ ^5+25^2 _|_ ^4+25^3/2 



<CL,ae'+''. (5.3) 



For the nonlinear part of the energy, Proposition 13. 31 gives 



{l~Cs'+'')l |a_r</ \^Pr<il + Ce^^-n I |a_| 



1+2(5n 



and the bound, |e ^'-^^ a^\^ = /^a 1^-1"' < Ci^r.Oi leads to 



2+2A- ^ ^1+25 



which is smaller than the error term for the linear part. This ends the proof of 



Remark 5.2. • The energy £-^{0.-) = £H,min + 0{tx/^) ■ Therefore the 

error given by (|5.2p is relevant, as compared with the energy scale of 

^ '^x£'T,Tnin: whcn 

e^^<Ciy,GTx, (5.4) 
with Ciy_G small enough, and accurate when 

e''<C,y,GrT. (5.5) 

Remember that the constants cg^ q o.i^d, Cg^ g depend also on 60, when 
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• It is interesting to notice that the worst term in the right-hand side of 
(|5.3p comes from the error of order 0{e^^'^^) in the Born-Oppenheimer 
approximation. There seems to be no way to get an additional factor r" 
with a > because the initial problem is rapidly oscillatory in the y- 
variable in a T^-dependent scale. This can be seen on the gain associated 
with the metric g^, for Tj, = 1 and > 0, which is simply (y^^p) or 
essentially 1 when p is small. 



5.2 Existence of a minimizer for Sr 



With the choice (|1.13p - (|1.14|) of the potential T4,t, the linear Hamiltonian HLi 
can be written 



Hl.u = -e''+''\^A + uoiq,T) 



f2^/l + T^x^ 
\ 

1/2 



with Wr{x, y) = 



For t < ^2p~' ^^'^ negative part ye 
Set 

HLin,t,+ = -e^+^^T^A + Uo{q,T) 



so that 



(l + r^a;2)2 1 + t^x^ ' 

part ^g2+2i5 ^(vj^-) is compactly supported. 



2vm^ 





uoiq,T)* 



\ 4 A' 



Gs.T 



2+25 



v{yjT^) t 



f2 



-2+25 



^). (5.6) 



Proposition 5.3. Assume e < £iv,G o,^d < t^^^q with t^^^q small enough. 
Then the infimum 

inf{£,(u), ueH\R'),\\u\\L2 = l} 
is achieved. Any element of Argmin £^ solves an Euler- Lagrange equation 

with the estimates 

|fe(V')l + |Av.|<Q,,Ge'+2^ , |l(l + r,|I?,nV2^,||^, <C,,,G, 

^ ^2+25 TxSH,min + ^'iv ,G^^^ +6^'') 
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Proof: From Proposition 15. 1[ we know that 

inf 4 (U) < e'+' V.£•r,n^^„ + Ci^ .Ge'+^*' 



II"IIl2=i 



For e < £iv,G and < ti^,-.Gi t smaller than ^^^a • Consider the decomposition 

(|5.6p for the energy £e{u) — f||u||^2 ■ By the same argument (convexity of the 
positive part and compactness of the negative part) as we used for £r in the 
proof of Proposition 14.61 (second step), a weak limit of an extracted sequence of 
minimizers in iJ^(M^) is a minimum for £^ on = 1} . 

A element ijj of Argmin satisfies 



G 



< {ip , HLin,t, + 1p) + 



2+2(5 



v{^.) t 



02 



2+2(5 



-Wr 



.2+2(5 



by recalling v > for the last line. This implies 



and by interpolation with ||u||£,6 < C||Vm 



2c; 



|2/3|i ||l/3 
L2 |lw||/2 



l^e = 0(t-^/'). The 



first inequality with ||V'||l2 = 1, gives + T^\Dq\'^)'^/'^ip\\L2 =0(1). 
The Euler-Lagrange equation 



implies 

|A| < 2 



(■0 , HLin,t, + i^) 



Ge.T 



2+2(5 



^'(^/T^•) t 



-2+2(5 



V') < c; 



II ^2+2(5 



Similarly, the lower bound £e{ip) > —2e^^^^ is due to Wr > —2. The upper 
bound oi £^{ijj) = £e.min comes from Proposition 15.11 and (|5.ip . □ 



5.3 Comparison of minimal energies between 8e and 

In this subsection, we specify a priori estimates for the minimizers of £e and 
compare the energies £e,min and £r,min without imposing relations between e^^ 
and Tx ■ This is not necessary at this level, if one uses carefully bootstrap 
arguments. 
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Proposition 5.4. Let V^^r be given by (|1.13p - (|1.14p and assume < t^^^g and 
£ ^ ^iv,G so that Ee admits a ground state according to Proposition\5^. The 
operator U is the unitary transform provided by Theorem \2.1\ and an element ip S 

Argmin £^ is written U 



'^^ I , with a = ( '^^ ] . Then, the estimates 



\\{l+T,\Dg\y^\\\L2<Cl,.^G, 



2+45 



2+45 



(5.7) 

(5.8) 

(5.9) 
(5.10) 



hold with right-hand sides which can be replaced by Ce^^G^^'^^^'Tx when e^^ < 



Proof: For d 



1 and a = 



the norms \\{l+T^\Dq\'^)'^/'^a\ 



L2 



and 11(1 + Tj.\Dq\'^y^'^d\\i^2 arc uniformly equivalent because 



Hence it suffices to estimate y^eDqU*ip . Remember that U = U{q,eDq,T,, 
with U € Su{l,gr;M2{C)), while G Su{{^p) , gr;^^) ■ With l|(l 
Ta;|£'qp)^/^'0||i2 < Ci^r^G in Proposition 15.31 simply compute: 

^eDgU*ij = eU*^Dgyj + [^eDg, U*^ i' = 0(e) in L''{R^; C^) , 

and divide by e for (|5.7p . 

In order to compare the energies, consider first the linear part by writing 



{'ip , Hunlp) - £ 



2+25, 



with 



BO — 



i?+ 



e 



■a+ 
~ a- 



-dl -{dy±i 



2V1 + T,X^ 



^ + 4^(l + r..Y/^ 
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Here it is convenient to write 



e TxtiBo ~ £ n ^ ' 



B = 



B+ 
S_ 



B± = B±{q,eDg,T,e) , 



(5.11) 



-V^)r + ^ . (5.12) 



B±{q,p,T,e) = t^pI + [y/T^Py ± T;i-rr-, , 

^ V J- + TxX 

With the notation x = x('''x£'^\Dq\'^) of Proposition 13.21 Lemma [5.51 below says 
in particular 

\\il-m\h<Ce,.,G [{a,m)+e'+''\\a\\l.] , 
\\^e\Dg\il-x)m,<Ci,.,G [{&,B&)+e'+^'Ml, 
Proposition 13.21 yields 



a, ( U*HMnU'e^^B 



p(l + T^x2)V2 



.2+4(5 I ^l+2S/~ 



a , Ba) 



For the nonlinear part of the energy, Proposition 13. 31 gives 

(l-Cei+2^)/ m^<[ \~at=( \a\^<{l + Ce'+'') [ 

JR2 JR2 JR2 Jr2 



with a = 



and a ~ 



coming from fe(V') = £'(£2+2'^) with e^+^-^VP^ = 0{e 



,2+25 Ti/ _/o^=-2+25n 



The bound /r^ IV'I'^ < QvCT"^ 



leads to 



2+2(5 



-^^1+26 



3+4(5 



which is smaller than the error term for the linear part. We have proved 



£,i^) - e^\a , Ba) - (5+ , 2(1 + TxX^ f^'~a+) - ^ \a\ 



e'+^^ +exe''{d, Ba) 



2S/~ 



With Seiij) = 8e,rmn (= 0{e'+'')), this giveS 



{1-Ce)e^^{h, Ba)+2{h+, {! + ■ 



„2^l/2- 



J \af<Ce^+^'+£, 



where all the terms of the left-hand side are now non negative. We deduce (|57 
and by bootstrapping 



\a\^<e^\a, B~a) + 



2+4(5 
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Additionally, / \a\'^ = J (|a_|2 + ja+p)^ > / \a^\'^ also gives 

< ||a- ||i2£r,)ri'in < Srici-) < Cg^^G ! 2+2S — ^'r.min + C'l^. q . 

With ||a_||i. = 1 - l|a+f = 1 + 0(e2+25) = l + 0(£!i) and £,,™„ = £^^,™„ + 
0{tx^^), this finally leads to 

^ ' X ' X 

□ 



Lemma 5.5. Assume B{q,p,T, 



T^|p|2+£r withr e Su{{^/^p),gT\M2{C)), 



take any x G C5^(R) and set B ~ B{q,eDq,T,e) and x — x{'^xA^q\) ■ Then 
there exists Eb^x ^ ^^'^ ^B,x > such that the estimates 



iu , Bu) + £ 



1+2.5 I 



-1+2(5 1 



11(1 - X)«lli^ + \\V^x\eD,\{l - x)u\\i. < Cb^x 
and ||(f - x)u\\l2 + \\^\eDq\{l - x)u\\l'^ < C'b^x I II^^IU" + ' ""II^^IIl^ 
hold uniformly w.r.t 6 £ (0,(5o] and e G (0,eB,x) • 

Proof: For xo G C^(R) sueh that xo ^ and xo = ^ around 0, the symbol 
Xo{txP'^) + S is an elliptic symbol in Su{{y/T^p)^ , gr ■ A^2(C)) . For e > small 
enough according to (xo, B), its quantization xo + B is invertible and its inverse 
belongs to OpSu{{\/T^p)^'^ , gT,.M2{C)) . Next we notice that in 



(i-x) = (i-x)' 



Xo + B 



Xo + B 



the last term belongs to OpNu,g^ if xo ^ X ■ AH the estimates are consequences 
of 



(1 - x) = {i-x)' 

withpe 5„(^^-iy^,.g,;X2(C)). 



Xo + B 



□ 



5.4 Adiabatic Euler-Lagrange equation 

As suggested by Proposition l4.71 or the last steps in the proof of Proposition l4.61 
an accurate comparison of minimizers requires some comparison of the Euler- 
Lagrange equations. We check here that the a_ component in Proposition 15.41 
solves approximately the Euler-Lagrange equations for minimizers of £r ■ Here 
the bootstrap argument is made in terms of operators instead of quadratic forms. 
In order to get reliable results, we now assume 

with ciy^G chosen small enough. 
With such an assumption we know: 
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• from (|5.10|) and \£r,min - £H,r,vin\ = 0{Tx^^), 

^-1^2+25 ^ c . / ^^2+25 

/ i ^ \ 

• When -0 g Argmin ^e^mm is written -0 = J7 y the L -norm 



of a is uniformly bounded, / \a\^ < Ct^^Ci according to ()5.8p . By apply- 
ing Lemma [3.41 this also gives / IV"!^ ^ C^v-.G- We also have Ija+Ula < 



n e.2+2(5_ 



• The Lagrange multiplier A.0, associated with -0 e Argmin £e,mim equals 
£e(0) + %^ / IV'I'' • Thus it is of order 0{e'^+'^^T^) . 

Proposition 5.6. Under the same assumptions as in Proposition \5.4\ and with 
the above condition e^^ < Ci^.^cTx, write a minimizer ip of Eg in the form ip = 

U f ^ I . Then the component a_ solves the equation 



H-a- + G \a-\'^a- ^ X,f,a- + re with \\re\\L2 < C^^,,G{£ -\ ), 

J Tx 

while ||a+i|L2 < Ct^, .ge'^+'^^t^ . 

The -norms of a and ip are uniformly hounded by Ci^-^g for p € [2,6] . 

/ I ^ \ 

Moreover if B is the operator defined in (j5.1ip - (|5.12p . a ~ y satis- 

fies 

\\B~a\\L^ <Ci^.^G£^T,,. (5.13) 

Remark 5.7. The relation of\\Ba\\]^2 with ||a||/f2 is given by 

[e^Tx\\Dl&\\ + e'WdU^] < \\m\\L^-+£^ML- < C [eV^p^^H + £2||a||^,] . 

Note that the L? remainder terms have a factor and not e^Tx ■ 

Proof: Playing with the Euler-Lagrange equation for -0, we shall first prove 
(|5.13p by using the same argument as we did for (a , Ba) in the variational 
proof of Proposition 15.41 and then use it in order to estimate \\r\\L2 . The Euler- 
Lagrange equation for i/' 

becomes 

Remember that Born-Oppenheimer Hamiltonian is given by 



-2-l-2<5^ rr _ ^2+25^ 

-- TxIdBO — £ Tx 



"H+e 
e~'^H-e' 



e^'B 
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by using the notations of (|5.1ip - (|5.12p . 

Let us consider first the hnear part after decomposing d into a = xa+(l— x)d, 
where the kinetic energy cut-off operator x = xi^'xl^Dgl'^) has been introduced 
in Proposition [ 



U*HL^nU~a = WHunUxa 

(I) 



+ [U* - U*)HunU{l - xYa + Lf^HuniU - Uo){l - x)d 



{II) 

(HI) 

The three terms are treated by reconsidering the computations done for Propo- 
sition By inserting a cut-off Xi i X ^ Xi i in XiU* HuntJxa with (1 — 
Xi)U*HLinUx € OpNu.g^, we get 



(/) = e'+''T.:,HBoXa + inL^(M^) . 

With e-i-25([> _ ^ OpSu[j^^^,gr;M2{'C)) and by using Lemma EH 
the second term is estimated by 

We write the third term as 

(///) = e''+^^T,HBo{l ~ xYa + e'+'^^T^VkU^ - x)~a + 2?p„t(f - x)~a , 

where Vkin and Vpot are defined by (|3.7p - p.8p . FoUowing the arguments given 
in the proof of Proposition l3.2l after these definitions, we get 

\\Vpot{l~x)~a\\L^ < Ce'+^'\\{l-x)a\\L-< 



< c 



e\\£^^I3h\\L2 +£ 



2+4(5 



- xYaU^ < C [£^+^^(1 - x)«IIl^ + e'+^'\\{e^AD,\)(l x)5|1l^] 



< C" 



e\\£^^I3h\\L2 +e- 



2+4(5 



Hence the Euler-Lagrange equation can be written 



with 



-2+2(5 



TxHboO- — ^ihO. — s 



2+2(5, 



-.,GIJ* (IVIV) +r 



< C 



e\\e^'Ba\\L^.+e^+'' 



From Proposition 15.31 with e^^ < ct^, we know HV-'IIl* < C . It can be trans- 
formed into ||a||i4 = |lC/*?/'||i4 < C with 



iiiv^r^'iiL^ <qiv^iii6 <c'iia|ii« 
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by applying Lemma 13. 4[ adapted for the metric gr like in the proof Proposi- 
tion [331 We start from the interpolation inequality 



\\f\\L<^<C\\D'j\\h\\f\\h< 



C 



I s 2 II 'X 



L4- 



~ 1 /O 1. 

By introducing the operator B = B{tx sq^ ri Dq^ r, e) with 



2 VvTT^ 

+ — 



Tx {Txe)Dy 



it becomes 



ll/IU« < 



l|s/III.II/lli4 + £S||/|II.||/|IE. 



— 1 — i 

The above relation with f ~ Tx ^ a{Tx ^ ■) which satisfies 

_i _i 
"llallLS , WfWh-^ = Tx "-\\a\\Li and \\B f\\L2 = \\Bh\\L2 



leads to 



< 



Cf, 



\\B^hTx^\~a\\l,+en~a\\l.,Tx^'\\~a\\l, 



{TxE) 

< Ci^,GTx^ \e-i\\Bd\\l, + l 



and finally 



s'/'rx\\e''B&U2+s'+-''Tx 



2+2(5^ 



With |A,0| = 0{e^+^^Tx), we obtain 

\\e^+^'TxHBoa\\L- < C \e\\e^' B~a\\L- + 

and recall 



2+45 I e.2+2<5_ 
\ I X 



_2+2AVic TT 
- J^BO 



e^^B+ + 2V1 + Txx^ 



According to Lemma [5781 below 

\\e^^B+~a+\\L2 + \\2^1 + Txx^h+\\L2 < C\\{e^^B+ + 2^1 + TxX^ya+\\L2 . 
We deduce 



2+25, 
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Plugging this result into the estimate of the remainder r gives 

Consider now more carefully the nonlinear term 

il* (IV'lV) = U*^ [|?7oa|2(J7oa)j + {U* - US) \\Uod\''{Uoii) 



(1) 



(2) 



u* 



\Ua\\U5.)-\Uo&\'iUoii) 



(3) 

By differentiating the relation J \ f\^ = J \Uof\'^ w.r.t /, the first term equals 
By semiclassical calculus in the metric g-r, the operator y/r^eDqUo equals 



y/r^eDgUo = Uoy/r^eDg 



TxeDq, Uo 



^ = U^V^eDq + 0{e) , 



where the remainder estimate holds in C^)) . We have already proved 

— and Lemma [3^ leads again to 



With £-i-25(^ _ (j^^ ^ OpSu[l,9r]M2[<C)), this gives 

11(2)11 <Cei+2^ 

For the third term, we use 



II \Udi\'{Uh)~\Uo~a\'{Uo~a) 1^2 < Co||(C/ - C/o)5||l« l|C/5|li6 + ||[/o5||i6 



< Ce 



1+2S 



Wc have proved 



which is even better than the estimate for the linear part. 
For the a_ component, we get 



H 



BO 



\ . ^/i- |2 , I- |2\~ _ -^V- 



-25 



+ G(|a_|^ + |5+n5_ 



and it remains to estimate the term |a+pa_ . The first line of the system may 
be written 
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Taking the scalar product with a+, with Ar = O {e'^'^'^^ t^) and > 0, gives 

■2S 

\~a+\' < C[\\a+\\h + {s + — )||a+|U=] < C'e'+^'r, . 
Then the same argument as in the estimate of || 511^6 gives 



l|5+IU« < C 
< C 



e-^||i3+a+||2.||5+||l4 + ||a+||l.||a+|||. 



2+25 „ 



-2+25, 



< C"Ti{e'+''Tx)^ <C" 
Again with ||a||L6 < C, we deduce 



4 + 4^ i 



The final result is just a transcription in terms of a . 
Lemma 5.8. Let he the symbol 



□ 



2Vl + r,a;^ 



/2 



introduced in (|5.12p . iJy setting _B+ = B^[q,eDq,Tx,e), the operator A = 
e'^^B+ + 2VTTx~^2 45 self-adjoint with D{A) = [u & L'^{B?), Au G L'^{R'^)} 
as soon as e < Eq, with Eq independent of . Moreover the inequality 



Vu e D{A), \\e''^B+u\\l2 + ||2v/l + r,a;2u||i2 < C||AM||i2 
ftoZds wit/i a constant C independent of (e, r^) . 
Proof: The operator A can be written 

A = ao{q, e'+^D„ Tx) + E^+^ai{q, e'+^ D,, r,) + E^+^^a2{q, r,) , 
with ak e 5„((V^p)2-'= + (^T^a;) 5r) and 

ao{q,p,Tx) + 2v/l + ri:a;2 . 

Therefore the operator A is elliptic in OpSuHy^p)'^ + {y^x) , gr) and the result 
about the domain follows with 



\\e'-'B+u\\l2 



for all u E D{A) . We conclude with 

2\\u\\l.<{u,Au)<\\uU4AuU. 



due to S+ > . 



□ 
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5.5 End of the proof of Theorem 11.21 

Assume G,lv > be fixed and S G (0,(5o] . Although we dropped Sq in our 
notations, all the constants in the previous inequalities depend on (G, iy , 6q) . 
We assume now t^; < Tgy^cSo, £ < ^iv.G.So and 

Proposition 15.41 says 

If . _ p2+25 c I < ^^2+45 < 2+25 1 

and we recall ^ 

\^T,min ^H,min\ ^ CTx 

by Proposition [ 



When ■(/; = [/ ' y , Proposition 15.61 savs 

\e a- I 

\\a+\\L^.<Ce^+^'T, , \\a\\H2<- and IjalU* + |la|Ua < C , 
while a_ solves the approximate Euler-Lagrange equation 

H_a^+G \a^\^a_=X^a-+rs with \\r^\\L2 < Gi^^ci^ ~l )• 

J Tx 

For u = ||a_j|^2^a_ the energy £t{u) satisfy 

\Et-{u) — (fr/miTil 1^ GTx 

and the above equation becomes 

H^u + G / \u\^u = \^,u + 0{e + — ) . 

./ Tx 



We conclude by referring to Proposition 14.71 applied to u and then renormal- 
izing for a_; For x (xi,X2) with xi e Cff (R^), x? + X2 = 1> Xi = 1 in a 
neighborhood of 

\\X2{tI .)a^\\L2 <Grl 

i 2 

\Sr{Xl{Ti a-) - £H,rnin)\ < Gt^ 

d-Hiixii^x .)a_, Argmin £h) < G{tx ^ + s) , vi^-^g S (0, -] . 

For the L°°-estimates of a+ and dicxi (xi(r^^ .)a_, Argmin Eh), we simply use 
the interpolation inequality 

IklU- <Ch||J.||A^i||J, <C||ii||J.||ii|||., 

valid for any u € i7^(]R^) (write u{0) = '^{Oj0)^^ cut the integral according 
to 1^1 ^ R, estimate both term by Cauchy-Schwartz with u = j|p-(|^pw) when 
ICI > Ri and then optimize w.r.t R). 
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6 Additional comments 



We briefly discuss and sketch how our analysis could be adapted to other prob- 
lems. No definite statement is given. Complete proofs require additional work, 
which may be done in the future. 



6.1 About the smallness condition of e w.r.t Tx 

In our main result, Theorem II. 2 [ condition (|1.18p is used, namely 



,25 < 



Cs ■ 

One may wonder whether such a condition is necessary in order to compare the 
minimization problems for £^ and e^^^^TxSh ■ When comparing the minimal 
energies in Proposition [531 we found 



|c . _ 2+25 . I ^ ^^2+45 

while we know that SH.min = £ h ,min{f-v ^ G) is a positive number independent 
of Tx and e . Hence it seems natural to say that e^^ ^ r^;, at least, is required 
to ensure that s^^'^^ Tx£ h ,min is a good approximation of £e,min ■ The error is 
made of three parts: 

• the error term for the Born-Oppenheimer approximation in the low-frequency 
range given in Theorem 12. Ij 



• the error term coming from the truncated high frequency part; 

• the non linear term. 

The non linear term is '•^'^ ^ "^"^ / IV"!^: so that a small error in ||'0|li4 will give 
a negligible term w.r.t e^^'^^Tx£H,min ■ The question is thus mainly about the 
linear problem. If one looks more carefully at the error term of Theorem l2.11 it 
is made of the term 

~ — E+~ E — -^^-^"'^ y^-^) 

according to Proposition 12.61 and of terms coming from the third order term 
of Moyal products. The function /e is in our ease fe{p) ~ £^^Txp'^"f{TxP^), 
P = (PxTPy), k,£ € {x,y} , and the factors Xx and Xy computed in the proof 
of Proposition 13.11 are at most of order . Hence the quantity (|6.ip is an 

0{e'^~^'^^Tx) which is again negligible w.r.t e^^^^Tx£H,mm ■ By considering the 
higher order terms in the Moyal product, the fast oscillating part of the symbol 
w.r.t y, at the frequency "^^i deteriorates the estimates: although there are 
compensations with the slow variations w.r.t py, always multiplied by y/r^, only 
an e'^ factor without Tx appears in the fe-th order term. 

Hence, computing the higher order terms, at least up to order 3, in the adiabatie 
approximation and then considering the question of the high-frequency trunca- 
tion, is a way to understand whether the smallness of e w.r.t Tx is necessary. 
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6.2 Anisotropic nonlinearity 

Our work assumes an isotropic nonlinearity. A more general nonlinear term 
would be 



ailV'il'' + 2ai2|?/'ini/'2p + Q;2|V'2|'' dxdy, ip 



Our case is ai = a2 = ai2 = 1. Let tp — Ucf), with the unitary transform 
U = Uq + e^+'^^R, (or conversely (j) = U*Tp). Then the same arguments as in 
Subsection 13.31 will lead to 



ail^il^ + 2ai2|Vi|'K'2|' + a2|^2|^ dxdy 

ai|^?|4 + 2ai2|^?n^'°l' + «2|V'2°I' dxdy^ (1 + 0{e'+^')) , 
after setting ~ '^a{<l)4'{Q) every q = {x,y) . In our case 

uo{x,y) 

with 6 = Oi^/r^x) , iy9 = . The point-wise identities 



C Se''^\ 



C = cos ( - I , S' = sin ( - ) , 



I^?P + I^2°P = I</'1P + |02P 

- 1^2 1' = cos(0) - |(/)2|') + 2sin(0) Me((/.ie»^</)2) . 



lead to 



I ;0|4 , o I ;0|2| ,0|2 I I ,0|4 + 2ai2 + 0-2 

ailV'il +2ai2|-0i| 1^2! +a2|V'2l = ] 



(I0l|' + I02|' 



Oil - "2 



(|0iP + |02p) cos(0)(|0i|2 _ |02p) + 2sin(0)Re(0ie^^(/)2) 



ai — 2ai2 + a2 



cos(0)(|0i|^ - |02r) + 2sin(0)Me(0ie'v02) 



"1 cos (-) + a2 sin*(-) + — sin2(6') 



I'/'il' 



+ 



aisin4(-) + «2COs4(-) + ^sin2(0) 
^^i±^sin2(0) + «i2(l + cos2(0)) 



i0in02 



+ [("1 - + (ai - 2ai2 + a2)cos(0)]sin(6')|(/ii|2Re(<^ie*'^02) 
+ [(ai - a2) - {ai - 2q;i2 + 02) cos(6')] sin(6')|(/)2p Ke(0le^(^2) 
+ (ai - 2ai2 + a2)sin^(6')Re ((he^(j>2) ■ 



At least three points have to be adapted from the previous analysis: 
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1) When we take a test function ip = U I -i-M— ] the energy Seii^) will be 



close to e^^^^T^ST-la^), with 



Er{a^) = (a_ , H^aJ) 
G 
+ 2 



ai cos'*(-) + a2 sin''(-) + — sin^(6') 



a_ 1^ dxdy , 



and cos(0) = /'f/^^ , • Hence before taking the limit — > we have a 
position dependent nonlinearity. This will induce another error term when 
comparing with the energy f/f(a_) = (a_ , iJ^y a^) + '^^"^"'""g^"'"^"^^^ /|a_|^ 
in the limit r^, — > . 

Another possibility consists in considering the case when |a2 — Q!i| + |ai2 — 
Qfil is small as (e,Ta;) — )■ (0,0) . The energy £^('0) j written as, 



G 

+ 2 



4 (q!i2 — Cki) . 

ai + (a2 - ai) sin^(2) + ^ ^ W 



will converge to = (a_ , Hi^.a^) + J |a_|^ . 

2) The existence of a minimizer for and the variational argument showing 

that ||a+||?2 = C'(e^+''*) + Se.min when tp — U i ^ ) is a ground 

state for 5^ will be essentially the same as in the isotropic case. 

3) The analysis and the use of the Euler-Lagrange equation for ground states 

of fg, like in Subsection 15.41 will certainly be more delicate because it will 
be a system, and the vanishing of the crossing terms have to be considered 
more carefully. 

6.3 Minimization for excited states 

One may consider like in |DGJO| the question of minimizing the energy, for 
states prepared according to ■0+ local eigenvector of the potential. Two things 
have to be modified in order to adapt the previous analysis: 

1) The space of states, on which the energy is minimized has to be specified. 
The unitary transform U introduced in Theorem 12.11 provides a simple 
way to formulate this minimization problem: set = Ran UP+ with 

P+ = \ \ I and consider 







inf feW. 

, l|V'll=i 
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2) In order to get asymptotically as e — >■ 0, the same scalar minimization prob- 
lems with Er and £h^ the external potential T^.t has to be changed. It 
must be now 



The analysis of this problem is essentially the same as for the complete minimiza- 
tion problem. It is even simpler because the unitary U is directly introduced. 
A slightly different question is about the minimization of the energy £^ in the 
space J-^ = RanC/o^+, but the accurate comparison between U and C/q widely 
used through this article would lead to similar results. 

Possibly this extension can even be generalized to higher rank matricial poten- 
tials with eigenvectors ipi, . . . jipN, for states modeled on any given ipk . 



6.4 Time dynamics of adiabatically prepared states 

Nonlinear adiabatic time evolution has been considered recently in jCaFe] . Note 
that our problem is slightly different because, we are considering a spatial adi- 
abatic problem, but some techniques may be related. 

When tpo ^ U ^ (resp. ipa ^ U (^^'^^ ), the question is whether the 
solution ip{t) to 



a-,oJ V 





remains close to U [ ^ ^{t)) (^"^^P" {^^Q^^^ ^'^^^ 

idta-=e^+^^Ta:H-a- + — -\a^fa- , a_ = 0) = a_,o , 

resp. idta^ = e^'^^^TxH+a^ H — ^|a+|^a+ , a+(t = 0) = a+.o . 

More precisely, the question is about the range of time where this approximation 
is valid: what is the size of w.r.t e such that sup^grj.^ rpj \\'ip(t)^U ' 



(resp. supjgf_j.^ ll'0(^) ^ ( '^^Q^'^ ) ID remains small. 



H+ 



Since the approximation of U*HLinU by - is good in the low fre- 

quency range, a natural assumption will be that the initial data are supported 
in the low frequency region 
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for some compactly supported x ■ Then the question is whether the norm of 
(1 - x(eV^D2))^(-^) remains small for t e [0,T^] for x G Cg"(R), x ^ Then 
the two last parts of Section [31 concerned with the high frequency part and the 
effect of U on the nonlinear term, have to be reconsidered. 

Note that the adiabatically prepared state with ipo = U ^^q^^ ^re probably 

not stable for very long time, = ©(et), because the characteristic set 

Cx = {{q,p) e R^ det {s^'t^p^ + V,Aq) + M{q) - A) = 0} 

contains two components when A > minii'+, one corresponding to the higher 
level of M{q) with |pp < C(A), and another one for the lower level of M{q) 
but with large p's. This means that a tunnel effect will occur between the two 
levels, so that adiabatically prepared states, with energies close to A > miniJ_|_, 
will not remain in this state for (very) large times. 



A Semiclassical calculus 

A.l Short review in the scalar case 

Consider a Hormander metric, g, that is a metric on R^"^, which satisfies the 
uncertainty principle, the slowness and temperance conditions (see |Horl [BoLej ) 
and consider ^-weights (slow and tempered for g) M, Mi, M2 ■ The symplectic 
form on M^J^ is denoted by a : 

d 

a{X,X') - ^p^g; - q,p" , X = {q,p),X' = {q',p') . 
i=i 

The dual metric 5'^ is given by g'^{T) = supy/_^Q ' g^[T') ^^"^ ^^'^ gain associ- 
ated with g is 

A(X) = mf (^^) ^ > 1 (uncertainty) . (A.l) 

In the scalar case, the space S{M,g) is then the subspace of C°° {M.p'^^;C) of 
functions such that 

yNeN,3CN sup \{Ti...TNa){X)\<CNM{X), 

gx{Ti)=i 

after identifying a vector field with a first-order differentiation operator. For 
a given N €N, the system of seminorms {pN,Ai,g)NeN defined by 

Pw,M,g (a) = sup sup M(X)-^\{Ti,...,TNa){X)\ 
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makes S{M, g) a Frechet space. 

For a G ) and £ > 0, the Weyl quantized operator a^{q,eDq) : 5(M'') 

iS'(IR'') is given by its kernel 

(q-q')-p ,q + q' . dp 



a'^{q,eD,){q,q') 



When no confusion is possible, we shall use the shortest notations 

a^{q,eDq) = a'^ = a{q,eDg) . 

When a € S{M,g), it sends iS(IR'^) (resp. S'{R'^)) into itself and the composition 
aioa2 makes sense for aj G S{Mj,g) . The Moyal product ai|l'a2 is then defined 
as the Weyl-symbol of aY{q,£Dq) o 02 [q.eDq): 



X\ — X^, — X 



>J-1 



X\ —X2 —X 



Xi — X2 — X 



+ e-^Rj{ai,a2,e){X) , (A.2) 



where Rj{.,.,e) is a uniformly continuous bilinear operator from S{AIi,g) x 
S{M2,g) into S{MiM2X~'\ g) (i.e. any scminorm of Rj{ai,a2,e) is uniformly 
controlled by some bilinear expression of a finite number of seminorms of ai and 
02). 

The three first terms of the previous expansion are given by 
aifa2 = 0102 + — [dp^aidqka2 - dqkaidp^a2] 

-J [(5p.,w«i)(^'^-,^«2) + id^qkya^)idl^p^a2) - 2(9p\.,.ai)(9,\.p,a2)] (A.3) 

+£^i?3(ai,a2,e) , 
by making use of the Einstein convention s^tj = sHj . 



Definition A.l. With the small parameter e € (0,£o), it is more convenient to 
consider the Frechet-space Su{M, g) of bounded functions from (0, £9) to S{M, g) 
endowed with the seminorms: 

PN.M.gia) ^ sup PN.M,g{a{£)) ■ 
ee(0,eo) 

The subscript „ stands for uniform seminorm estimates. 

The space of e- quantized family of symbols will be denoted by OpSu{M,g): 

(a(£) e Su{M,g)) ^ {a^'{q,eDq,e) e OpSu{M,g)) . 
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We shall give a variation of this definition in Subsection IA.4I below for pa- 
rameter dependent metrics. 

We recall the Beals-critcrion proved in jBoCh] (see |NaNi| for the e-dcpendcnt 

version) for diagonal Hormander metrics of the form g = ^p^lxy^ ~^ ^'^T^Y^ ' 

Set S) = {Dqi , . . . , Dyd , , . . . , q'^) and for E E N^** introduce the multi-commutator 
^•^s = nj='i acting on the continuous operators from 5(]R'') to S^W^) and 
the weight Me{X) = iip{X), ^(X))^ . 

Then the Beals criterion says that an operator A : iS(M'^) — belongs to 
OpS{l,g) if and only if 

ad| A e C{L^{R'^), H%Me)) 

for all E E , when the Sobolcv space H^{Me) is given by 

\W\\h-{Me) = UlE{q,eDq)u\\L2 . 

Moreover the family of seminorms {QN)NeN, defined by, 

Qn{A)= sup max e~^\\ad§,A{e)\\c{LHm'').H'{ME.g)) (^-4) 
ee(0,£o) \^\=" 

on OpSu{M, g), is uniformly equivalent to the family (P;v(a))Af6N on Su{M,g) 
after the identification A{e) ~ {q, eDg, e) . "Uniformly" means here that the 
comparison of the two topologies is expressed with constants independent of 
e e (0,eo)- 

Below is an example of a metric which satisfies all the assumptions and which 
is used in our computations 

with the gain function X{q,p) = (p)^ . 

It is convenient to introduce the class of negligible symbols and operators. 
Definition A. 2. An element a E Su{^,g), belongs to Mu.g if 

yN,N' E N, 3Cn,n' > 0, e-^a(e) E Su{X~^\g). 
Similarly, OpMu.g denotes the e-quantized version: 

{a^{q,eDg,e) E OpNu^g) ^ {a E Mu,g) ■ 

Combined with the following relation between cut-off functions it provides 
easy estimates from phase-space localization. 

Definition A. 3. For two cut-off functions XI7X2 & C^(R^''), < xi,2 < 1, the 
notation xi ^ X2 means that X2 ^ ^ in a neighborhood of supp xi ■ 

For example the pscudodiffercntial calculus leads to 

(Xi^X2)^(Vae5„(Af,5), (1 - X2)f af Xi e AA„,,) , (A.6) 
when the weig ht M satisfies M < CX^ for some C > . 
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A. 2 Applications to matricial operators 



Operator valued pseudodifFerential calculus has been studied in jBak] . When 
f} is a Hilbert space it suffices to tensorize the previous calculus with £(fl), 
which corresponds to the componentwise definition in A^„(C) when () = C" . 
The corresponding class of symbols associated with a Hormander metric g and a 
g-tempered weight M is denoted by S{M, g; C{t))) while the set of bounded fam- 
ilies in S'(M, 5; £(())) parametrized by e € (0,eo) is denoted by S'„(M, £(f))) . 
The asymptotic expansions (jA.2p - (jA.3|) of the Moyal product clearly holds (see 
|PST| for a presentation without specifying the remainder terms) if one takes 
care of the order of the symbols. For example, the two first terms of the expan- 
sion of a commutator [ai{q, eDq), a2{qT eDq)] are 

aif 02 - a2fai = [ai{q,p), a2{q,p)] 

£ £ 

+ -^{dpaidqa2 - dqaidpa2) ~ —{dpa2dqai - dqa2dpai) 

-f [R2{ai,a2,£) - i?2(a2, ai, e)] (A.7) 

with no simpler expression when the matricial symbols do not commute. 

When the Hormander metric q has the form ^ — h ^ the uniform Beals 

criterion also holds: In the scminorms Q{a) defined in (|A.4I) simply replace 
Li^iW^) by the Hilbert tensor product L'^{W^) ® t) and consider the operators D 
and Msiq, £Dq) as the diagonal ones D ® Id|, and Me{, £Dq) ® Id(, . 
Finally the Definition IA.2I of negligible symbols also makes sense for matricial 
symbols after replacing OpS{X~^ ,g) by OpS{X~^ ,g]C{i})) . 



We end this section with standard applications of the Beals criterion. 
Proposition A. 4. Consider a diagonal Hormander metric g — if-fx)'^ 

dp^ 

^.(^xy^ '^^'^ constant metric go = dq^ + dp^ . Assume that any f chosen in 
{ipj^ipj, j € {1, . . . ,d}} is a go-slow and -tempered weight such that f{q,£DqY 
belongs OpSuif, go) for any s e N . El 

Assume that A G OpSu{l, g', C{t))) is a family of invertible operators inC{L'^{W^)^ 
[)) and such that ||A^"'^(e)|| is uniformly hounded in w.r.t £ e (0,eo 
belongs to Op5ti(l, g; £(f))) . 

Proof: We start from the relation 

E CE,E'A-^^'^[{^^4^ A\A-^ 



ad| A-^ 



\E'\ = \E\,E'en\'^\ 



i=i 



Then A 



(A.8 



which holds in C{L'^{W^) ®\)) for all E e N^'^ . Hence the Beals criterion in the 
metric go says that belongs to OpS'„(l, £(f))) . In particular A'^^ belongs 
to C{H'^ {Me" )) for any iJ" e N'^'^ . This allows to apply the Beals criterion in 
the metric g and yields the result. □ 

■^This last condition is redundant after a possible modification of ipj and ipj if one refers 
to IBoChl . but easier to check directly in our examples than giving the general proof. 
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A. 3 PseudodifFerential projections 

An application of the Beals criterion says that a true pseudodifferential projec- 
tion can be made from an approximate one at the principal symbol level. This 
holds for matricial symbols. 

Proposition A. 5. Consider a diagonal Hormander metric g — ^\xY ^~ 

^.{XY with the same properties as in Proposition \A.4\ Assume that the operator 
ns OpSuil, g; satisfies 

(n o n - n) = e^R^ + e^k^ , 

with Rf, e Su{M,g; £{[))), R^ G Su{N, g; C{i))) v > > Q and M,N < I . 
Assume additionally that there exist Xi x' G Su{'^,g] < X ^ x' !i 1 such 

that X ^ x' x'Rfj. = . 
Then for ei < Sq small enough, the operator 

is well defined for £ € (0,£i) and satisfies 

PoP^P m C{L^{R'^)®t}) , 
and (^P-tlj = ef'tlf, + £"11^ . 

with e Su{M,g;C{t))), e Su{N,g;C{t))) and o x € OpMu,g ■ 

Proof: The first result concerned with the definition of P is a direct ap- 
plication of the simple general result in Lemma IA.6I applied with T ~ II and 
"H = L'^{M.'^) (g) f) . Note that our assumptions v > 1 and M < 1 ensure that 
||n^ — n|| < i as soon as £ < £i with £i small enough. 
Writing 

fl-P=(n^ -Ii){Ai-Ao) 

with 

A.^^ I {z-iiy' dz, {0,1}, 

reduces the problem to proving 

A,eOpSu{l,g]Cm or (z - fl)-i e Op5„(l,g;£(t))) , 

when {\z — a;| = 1/2} with uniform bounds. But this was proved in Proposi- 
tion IA.4I as a consequence of the Beals criterion. This ends the proof. □ 
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Lemma A. 6. Assume that in an Hilbert space %, the operator T G 
satisfies ||r^ — r|| < (5 < 1/4 and \\T\\ < C. Then there exists cs < ^ such that 

a{T) C {z e C,|z(2- 1)1 < (5} C {z G C,|z| < cs} U {z G CJz- 1| < cs} , 
niax|||(z-r)-i||,|z-l| = U <2^^. 



Moreover, the operator 

P = 

2i7r 



\z~l\ = l/2 



differs from T according to 

T-P={T^- T){Ai - Aq) = {Ai - A^){T^ - T) , (A.9) 

with Ai^^l {T ^ z)-\l - z)-^ dz (A.IO) 

2«7r7|3_i|^i 

and Ao^^l {T-zy^z^^dz. (A.ll) 



2i7r 



Proof: If z G a{T) then z{z - 1) G a{T{T - 1)) C {z G C, |z(z - 1)| < i} 
(Remember that |z(z — 1)| = |; means \Z —j\ = j with Z ~ {z~ ^Y)- Consider 
z G C such that |z — 1| = i, then the relation 

(r-z)(r-(i-z)) = z(i-z) + (r2-T), 

with |z(l - z)| > i and IjT^ - T\\ < 5 < \, imphes 

\\iT-zr^\\<\\T-{l-z)\\\\[z{l-z)+T^-T]-^\\<^^ for|z-l| = i. 

The symmetry with respect to z = i due to (1 - T)(l - T) - (1 - T) - T 

imphes also 

ll(T-^)-'ll<^ forlzH^. 

Compute 

T~P ^ ^1 [T{z-l)-^-{z-T)-'] dz 

n2 



(T- 1)P+ (T^ -T)^i 



with Ai ^ — I (T- z)-i(l - z)-i dz. 



In particular this implies to r(l — P) = (T^ — T)Ai while replacing T with 
(1 - T) and P with 1 - P leads to 

P-T = -r(l-P) + (T2-r)Ao with Ao^^l {T - z)-^ z-^ dz . 

2nT J\z\=i 
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Summing the two previous identities yields the result. 



□ 



A. 4 Extension to parameter dependent metrics 

Additionally to the scmiclassical (or adiabatic) parameter, we need other pa- 
rameters T = {t' ,t") € (0, 1]^ on which the metric g = gr depends. In general 
consider r G T C and a family of Hormandcr metrics (,9T)reT defined on 

Definition A. 7. The family of metrics {gT)Ti=T •said admissible if the uncer- 
tainty principle (|A.ip is satisfied and if the slowness and temperance constants 
Ci , 6*2 , N2 involved in 



can be chosen uniformly w.r.t t ^ T ■ 

Accordingly a family of weights {Mr)TeT "^^^^ be admissible if the slowness and 
temperance constants of Mr w.r.t g-r can be chosen uniform w.r.t t Cz T . 

The important point is that all the estimates of the Weyl-Hormander pscu- 
dodiffcrcntial calculus (see }Hor( iBoLe] ). including the equivalence of norms in 
the Beals criterion of [BoChj . occur with constants which are determined by 
the dimension d, the uncertainty lower bound (which is 1 here), the slowness 
and temperance constants. Hence all the pseudodifFerential and scmiclassical 
estimates, (operator norms or seminorms of remainder terms) are uniform w.r.t 
to (e,T) as long as the symbols, afc(£,T), k = 1,2, have uniformly controlled 
seminorm in S{Mr,k, 9t), w.r.t (£,t) G (0,£o] x T. 

The definition of symbol classes Su{M,g) with uniform control of seminorms 
w.r.t e € (0,eo] can be extended to admissible families {Mr, gT)T<£T ■ 

Definition A. 8. For an admissible family {Mr,gT)T£Tj set of parameter 
dependent symbol a{X,e,T), X G M.'^'^ , (e,r) G (0,eo] xt with uniform estimates 

V7V G N, 3Cn > 0, V(e, r) G (0, eo] x T , PN,AU,a. r)) < Cn , 

is denoted by Su{Mr,gT,^{^)), t £T ■ 

Equivalently the set of semiclassically quantized operators a{q,eDq,e,T) when 
the symbol a belongs to Su{Mr, gr] ^ih)) denoted by OpSu{Mr, gr', ^{h)) ■ 
The set of negligible symbols and operators associated with {gr)TeT with uniform 
estimates in Definition \A.S\ w.r.t t £T is denoted by Afu,gr '^'^^ OpNu^g^ ■ 

Proposition A. 9. For t — {t',t") G (0, 1]^ the family {gT)Te{o,i]^ defined on 



(temperance) 



(slowness) 




T'dq'^ 



+ T"dq"^ + 



T'T"dp^ 



9t = 




(V^p) 



2 
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is admissible. 

Proof: It is easier to consider the symplectically equivalent metric (use the 
transform {q' , q" ,p' ,p") (t'^ q' ,t"^ q" ,t'~^ p' ,t"~'^ p")) 

dq'^ , „ T'^dp'^ + T"^dp"^ 

'^^w^'' + — m — 

after setting (p)^ = 1 + t''^p''^ + t"^p"'^ . Firstly remember that the metric 

dq'^ dp2 
3(1,1) = ^ + ^'? 

is a Hormander metric. The metric is given by 

= ^dg'^ + + {q^dp'^ + dp"^ . 

Hence the uncertainty principle (jA.l[) is satisfied with 

A.(g',g",p',p") =min|(g')^,^| > (p). > 1. 

In order to check the uniform slowness and temperance of g-r, introduce the 
new variables Xr — {q' ,q" ,t'p' ,t"p") when X — {q' ,q" ,p' ,p") with a similar 
definition for Yr and . 



Slowness: Write 



(^gr,x{X _ y) < i-^ ^ - Y^) < ^) 



When Ci is slowness constant of 5(i,i), this implies 

±1 



which is nothing but 

\ ±1 

gT,x 

9r,Y 

Temperance: Write 

±1 ^\ ±1 



\9t,y[T)J V5(i,i),k,{Tt) / \ ' )' - 

when 6*2 and sre the temperance constants for • The problem is 
reduced to showing 

The expression of g^ gives with X = {q' , q",p',p") and T = {9' , 9", tt', tt") 

5rM)x.(r.) = {P)l9'' + {P)l6"^ + {q')W + r'"^'" < rr.xiT) 
owing to maxjr', r"} < 1 . 



73 



□ 



Acknowledgements: The authors are very grateful to Jean Dahbard for ex- 
plaining the mathematical questions arising from a two level atom in a light 
beam and for many detailed discussions that took place for the duration of 
this work. They also wish to thank G. Panati, C. Lebris and S.V. Ngoc for 
their occasional help. This work was initiated while the second author had a 
CNRS-sabbatical semester at CMAP in Ecolc Poly technique. They acknowledge 
support from the French ministry Grant ANR-BLAN-0238, VoLQuan. 



References 

[Aft] A. Aftalion. Vortices in Bose-Einstein Condensates. Progress in Nonlinear 
Differential Equations and their Applications, Vol 67 Birhkaiiscr (2006). 

[AB] A. Aftalion, X. Blanc. Reduced energy functionals for a three dimensional 
fast rotating Bose Einstein condensates. Ann. Inst. H. Poincare Anal. Non 
Lineaire 25 no. 2 (2008), pp 339-355. 

[ABN] A. Aftalion, X. Blanc, F. Nier. Lowest Landau level functional and 
Bargmann spaces for Bose-Einstein condensates. J. Funct. Anal. 241 no. 2 
(2006), pp 661-702. 

[AJR] A. Aftalion, R. L. Jerrard, J. Royo-Letelier. Non existence of vortices in 
the small density region of a condensate Jour. Funct. Anal., vol. 260, no. 
8 (2011), pp 2387-2406. 

[Bak] A. Balazard-Konlein. Calcul fonctionnel pour des operateurs /i-admissible 
a symbole operateur et applications. ph-D, Universite de Nantes, (1985). 

[BCR] J. Bochnak, M. Coste, M.F. Roy. Geometrie algebrique reelle. Ergebnisse 
der Mathematik und ihrer Grenzgebiete (3), 12 Springer- Verlag (1987). 

[BDLM] J. Bolte, A. Daniilidis, O. Ley, L. Mazet. Characterizations of Lo- 
jasiewicz inequalities and applications: subgradient flows, talweg, convex- 
ity. Math. Oper. Res. 36 no. 1 (2011), pp 55-70. 

[BoCh] J.M. Bony, J.Y. Chemin. Espaces fonctionnels associes an calcul de 
Wcyl-Hdrmander. BuU. Soc. Math. France, 122 no 1 (1994), pp 77-118. 

[BoFo] M. Born, V.Fock. Bcweis dcs Adiabatensatzcs. Zcitschrift fiir Physik 51 
(1928), pp 165-169. 

[BoLe] J.M. Bony, N. Lerncr. Quantification asymptotique et microlocalisation 
d'ordre superieur L Ann. Scient. Ec. Norm. Sup., 4*^ serie 22 (1989), pp 
377-433. 

[BoOp] M. Born, R. Oppenheimer. Zur Quantentheorie der Molekeln. Ann. 
Phys. (Leipzig) 84 (1927), pp 457-484. 



74 



[BrOs] H. Brezis, L. Oswald. Remarks on sublinear elliptic equations, Nonlinear 
Analysis 10 (1986), pp 55-64. 

[CaFc] R. Carles, C. Fermanian. A nonlinear adiabatic theorem for coherent 
states. Nonlinearity 24 no. 8 (2011), pp 2143-2164. 

[CFKS] H.L. Cycon, R.G. Froese, W. Kirsch, B. Simon. Schrodinger operators 
with application to quantum mechanics and global geometry. Texts and 
Monographs in Physics. Springer Study Edition. Springer- Verlag (1987). 

[Chi] R. Chill. On the Lojasiewicz-Simon gradient inequality. J. Funct. Anal. 
201 no. 10 (2003), pp 572-601. 

[Coo] N.R. Cooper. Rapidly Rotating Atomic Gases Advances in Physics 57 
(2008), p 539. 

[DGJO] J. Dalibard, F. Gerbier, G. Juzeliunas, P. Ohberg. Artificial gauge 
potentials for neutral atoms. Rev. Mod. Phys. 83 (2011), p 1523. 

[Fet] A.L. Fetter. Rotating trapped Bose-Einstcin condensates Rev. Mod. Phys. 
81 (2009), p 647. 

[LCGPS] Y.J. Lin, R.L. Compton, K.J. Garcia, J.V. Porto, LB. Spielman. Syn- 
thetic magnetic fields for ultracold neutral atoms Nature 462 (2009), p 628 

[GCYRD] K.J. Giinter, M. Cheneau, T. Yefsah, S.P. Rath, J. Dalibard. Prac- 
tical scheme for a light-induced gauge field in an atomic Bose gas. Phys. 
Rev. A 79 (2009), p 011604(R). 

[HaJe] A. Haraux, M.A. Jendoubi. The Lojasiewicz gradient inequality in the 
infinite-dimensional Hilbert space framework. J. Funct. Anal. 260 no. 9 
(2011), pp 2826-2842. 

[HiPr] M. Hitrik, K. Pravda-Starov. Spectra and semigroup smoothing for non- 
elliptic quadratic operators. Mathematische Annalen, 344 (2009), no. 4, pp 
801-846. 

[Hor] L. Hormander. The analysis of linear partial differential operators. 
Springer Verlag (1985). 

[Hua] S.Z. Huang. Gradient inequalities. With applications to asymptotic be- 
havior and stability of gradient-like systems. Mathematical Surveys and 
Monographs, 126. Am. Math. Soc. (2006). 

[Kat] T. Kato. Perturbation theory for linear operators. Reprint of the 1980 
edition. Classics in Mathematics. Springer- Verlag, (1995). 

[LiSe] E.H. Lieb, R. Seiringcr. Derivation of the Gross-Pitaevskii Equation for 
Rotating Bose Gases. Commun. Math. Phys. 264 (2006), pp 505-537. 



75 



[Loj] S. Lojasiewicz. Une propriete topologique des sous-ensembles analytiques 
reels. Les Equations aux Derivees Partielles, pp. 87-89, Editions du centre 
National de la Recherche Scientifique (1963). 

[MCWDl] K.W. Madison, F. Chevy, W. Wohlleben, J. Dalibard. Vortex for- 
mation in a stirred Bose-Einstein condensate. Phys. Rev. Lett. 84 (2000), 
p 806. 

[MCWD2] K. W. Madison, F. Chevy, W. Wohlleben, J. Dalibard. Vortices in 
a stirred Bose-Einstein condensate Jour. Mod. Optics 47 (2000), p 2715. 

[MaSo] A. Martinez, V. Sordoni. A general reduction scheme for the time- 
dependent Born-Oppenheimer approximation. C.R. Acad. Sci. Paris Ser. 
I, 334 (2002), pp 185-188. 

[MaSo2] A. Martinez, V. Sordoni. Twisted Pseudodifferential Calculus and 
Application to the Quantum Evolution of Molecules. Mem. Amer. Math. 
Soc. 200 no. 936 (2009). 

[NaNi] F. Nataf, F. Nier. Convergence of domain decomposition methods via 
semi-classical calculus. Comm. Partial Differential Equations 23 no. 5-6 
(1998), pp 1007-1059. 

[NeSo] G. Nenciu, V. Sordoni. Semiclassical limit for multistate Klein-Gordon 
systems: almost invariant subspaces, and scattering theory. J. Math. Phys. 
45 no. 9 (2004), pp 3676-3696. 

[Nie] F. Nier. A propos des fonctions theta et des reseaux d'Abrikosov. 
Seminaire Equations aux Derivees Partielles. Ecole Poly technique, 2006- 
2007, Exp. No. XII. 

[PST] G. Panati, H. Spohn, S. Teufel. Space adiabatic perturbation theory. 
Adv. Theor. Math. Phys. 7 no. 1 (2003), pp 145-204. 

[PST2] G. Panati, H. Spohn, S. Teufel. The time-dependent Born-Oppenheimer 
approximation. M2AN 45 no. 2 (2007), pp 297-314. 

[Sim] L. Simon. Asymptotics for a class of non-linear evolution equations, with 
applications to geometric problems. Ann. of Math. 118 (1983), pp 525-571. 

[Sjo] J. Sjostrand. Parametrices for pseudodifferential operators with multiple 
characteristics. Ark. fiir Mat. 12 (1974) pp 85-130. 

[Sor] V. Sordoni. Reduction scheme for semiclassical operator-valued 
Schrodinger type equation and application to scattering. Comm. Partial 
Differential Equations 28 no. 7-8 (2003), pp 1221-1236. 

[Tay] M. Taylor. Palatial Differential Equations III, Nonlinear Equations. Ap- 
plied Mathematical Sciences Vol. 117, Springer (1997). 



76 



